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I. INTRODUCTION 


The concept of observability introduced by Kalman in 1960 [Ref. 1], is, along 
with the dual concept of controllability, one of the most significant ideas to come out 
of the state space approach to control theorv. The concept of observabilitv presented 
in most of the literature available today, gives only an answer to whether the svstem 1s 
observable or not. However. another question was raised by Brown in 1966 [Ref. 2). 
opening the issue of how observable a given system is. Are all the states equally 
observable? If not, what are the relative degrees of observability of the state variables? 
Since 1966, this issue has been investigated in [Ref. 3]. [Ref. 4]. [Ref. 5] and [Ref. 6] 
among others. In the majority of theses investigations the goal has been either the 
state-space reduction of the model or frequency domain model reduction. Silverman 
has probably been one of the most active researchers on this subject, using, in several 
of his papers, the technique of rotating the coordinate system in such a way that a so 
called “balanced realization” is obtained. By studving this new realization the “nearly 
redundant” states can be distinguished from the most important ones. providing the 
required information to reduce theierdemo! mie moder 

In this thesis the problem of how observable the states are is again studied using 
techniques different from the ones mentioned above. Two graphical methods are 
developed: the first one 1s relatively simple, but rather helpful in understanding the 
concepts involved in ranking the state variables according to their relative degrees of 
observability. The second graphical method, based on some ideas given in [Ref. 7], is 
more elaborate. Both graphical methods are compared against a third method 
described in [Ref. 7] whose development is explained in this thesis. 

The tests to find the relative observabilitv for each state variable provide the 
information to develop a method that enables the improvement of the degree of 
observability of the least observable state. Although all the basic derivations assume a 
single output regulator problem. the generalization to multi-output non-regulator 
svstems is also made. This thesis is, therefore, divided as follows: in Chapter II a brief 
and general definition of observability is given followed by two single output numerical 
examples. Chapter III begins with a definition of relative observability followed by the 


presentation of methods to evaluate the relative degree of observabilitvy of a system. 


Three algorithms are presented leading to a Fortran computer program whose listing 1s 
in Appendix A. In Chapter IV an algorithm to improve the relative degree of 
observability of a weakly observable state is developed for the single output regulator 
problem. This algorithm leads to a Fortran program that is contained in Appendix B. 
The generalization for multi-output servo systems is also discussed. Finally in Chapter 
V the performances of the three methods of Chapter III are compared as well as the 
performance of the algorithm developed in Chapter IV. Some particular examples used 


with the computer programs are presented and analized in some detail. 


If. THE CONCEPT OF OBSERVABILITY 


A. DEFINITION 

A system is said to be completely observable if, from the measurements of the 
output v(t) over a finite period of time, t €(0.ty]. every initial state can be determined 
(Ref. 8: p. 764]. [Ref. 9: p. 107]. In other words, if every transition of the states afiects 


every element of the output vector, the system 1s completely observable [Ref 8: p. 764]. 


B. THE STATE-OBSERVABILITY PROBLEM 


Consider a linear time-invariant svstem described by equations 2.1 and 2.2 


X(t) = Ax(t) + Bu(t) (oun 


¥(t) = Cx(t) (2.2) 


where: 

x(t) is a nx! column vector of the state-variables. 

x(t) is a nx!l column vector of the time derivatives of the state variables 
A IS a nxn matrix. 

B is a nxp matrix 

uis a pxl column vector of the inputs to the system. 

\(t) 1s a mx1 column vector of the output variables of the system 


Ci8 4 mn en 


The effect of the input signal u(t) can always be determined, if x(0) can be 
computed from the measurements of the output. 1e., if the system is completely 
observable. Therefore, there 1s no loss of generality by assuming that u(t)=0 [Ref. 9: p. 
107]. Under this assumption, to evaluate a svstem’s observability a matrix Q, called 


the observability matrix, must be formed. [Ref. 10: p. 108]. 


Each of the submatrices separated by spaces in equation 2.3 is nxt. Therefore. 
there will be n columns in the Q matrix; this matrix will be square (nxn) when there is 
onlv a single output. [f the n columns of Q are all linearly independent. 1.e., if the nxn 
matrix Q is nonsingular the initial state vector x(0) can always be found from the 
measurements of the output [Ref. 10: p. 105]. On the other hand, if Q is singular. (not 
of full rank). some columns of Q will be linearly dependent on others, and only special 


choices of v(t) will produce the required solution [Ref. 10: p. 105]. 


Cab NANER LE 
1. Observable system 


Consider the single output linear time invariant system described by: 


x(t} = Ax(t) + Bu(t) 


eS ext) 
Where: 
1 2 0 0 C= E 0 1 | 
A=]2 3 0 B= 10 
Ory Lr l 


The matrix Q 1s equal to: 


eegs 
O= 0 ele 
Lo 


This matrix Q is nonsingular, ie.. of full rank, therefore the system 1S 
observable. 
2. Non-observable system 
Consider, again the same state representation where the matrices A, B and C 


are as follows: 
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The matrix Q 1s equal to: 


O- 0020 
Q=]0 22 
ees 


In this case the matrix Q is singular, 1.e., its determinant is equal to zero. Its 
rank is also less than three, (actually two), because the maxinium submatrix of Q with 
deternunant different from zero is a 2x2 matrix. Therefore, in this case there is one 
state that 1s not observable. The identification of that non observable state is not 
possible just by performing the previously described test, where the answer “ves-no” is 
given, With no indication of how observable the system 1s and which are the most 
observable states. In the following Chapters methods will be discussed and presented 


where answers to these questions will be given. 


Il. THE CONCEPT OF RELATIVE OBSERVABILITY. 


A. DEFINITION 

The problem of relative observability deals with finding an algorithm or figure of 
merit for each state-variable, which will reflect how observable the state-variable is. 
Measuring y(t) at times t = ty. ty,...t,, One wants to be able to evaluate which state 
components can be most accurately determined from the n set of measurements of the 
output. In other words, how an error in the measurements is reflected in determining 
the initial states? If a small error in the measurements produces a large error in the 
calculated state, this state 1s. therefore, difficult to determine from the measurement of 
the output. The output has a “fuzzy image” of the state that 1s “hidden” by the large 
error. If a small error in the measurements produces a small error in the evaluation of 
a particular state, this state will be evaluated more accurately than the others from the 
measurements of the output. Due to the small error the “image” of the state at the 


output is less “fuzzy”. 


B. METHODS FOR EVALUATING RELATIVE OBSERVABILITY 

The three methods that will be described for evaluating the relative degree of 
observability were developed based on the concept given previously: the projection of a 
measurement error in computing the initial state. 

Starting with equations 2.1 and 2.2 and assuming the input to be zero. equation 
2.2 is differentiated n-!| times and equation 2.1 1s substituted as shown below 
Per t0mpe 106]. 


¥(to) = CX(tg) 3.1) 


Egmont ty) CAX(tH) 
itp) GAK(to) = CA-x(tp) 


lty) = CA! x(tp) 


Is 


The result is n equations with n unknowns. It can be seen that the coeflicient 
VEClONS a Gn © an CA’... CAM) are row vectors with the same elements as the 
respective column vectors cl. Alclyalyecl ualyalel of the QO Mmarins. Wiis. 
the elements of the columns of the Q matrix are the coefficients of the equations 


relating the initial states with the measured quantities. Equation 3.1 can be written as 
-oly ) 14 
V(t) = Q°X(ty, (3.2) 


Where vy(tg) is the column vector of the output and its n-1 derivatives evaluated at t = 
ty: 

The fact that the output equation 2.2 must be differentiated n-1 times to obtain n 
equations and n unknowns, implies that the output has to be measured n times. If the 
system has only one output, as in this case, the state vector, x(tg), can be found by 
inverting the ql Matric sbheinversesen on onlv exists if its determinant is nonzero or, 
in other words, if the rank of the matrix is equal to its order. However, if the 
coefficients of any pair of equations are approximately proportional, the determinant 
of ql will be a small number and, from a practical point of view, it will be difficult to 
invert the matrix and solve the system of equations. A matrix that has a small number 
for its determinant 1s an “ill-conditioned” matrix. Since a zero value for the determinant 
occurs only for a singular matrix, it can be said that an ill-conditioned matrix is nearly 
singular. 

The solutions of ill conditioned systems are very sensitive to small changes in 
either the coefficients of the matrix or the values of the forcing functions. Therefore, if 
the rows of gl are nearly proportional, a small error in in the measurement of yq(tg). 


will produce a large error in x(ty) due to the high sensitivity of ill-conditioned systems. 


C. UPPER-BOUND ERROR 
|. Theoretical background 
The degree of independency of the rows of qi can be used to quantifv the 
degree of observability for each state variable [Ref: 7: p. 12]. If two or more of the 
rows of Ql are almost linearly dependent the entries of (Q hy} will be large numbers. 
In defining the degree of observability for each state variable using the upper-bound 
error, the reciprocal of the increase in the error of the calculated state-variable over the 


observation error will be used [Ref. 7: p. 17]. The reciprocal will be used because a 
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small number associated with a low degree of observabilitv is more intuitively 
appealing than a large number. 

Referring to matrix equation 3.2, normalization of the rows of on will be 
performed. The goal of this operation is to obtain a maximum value for the degree of 
observability equal to one. To accomplish the normalization the elements of q(t) are 
divided bv the length of the corresponding row vector of the Qi matrix. The 
normalized vector is defined as ¥n(to) and the normalized on as P. Equation 3.2 can 


now be written as: 


¥,{to) a Px( to) ( 


(oJ 
Us 


The vector y,(tg) consists of the actual value of the measurements v,(ty) and 


an error e.,,. thus, equation 3.3 can be written as : 
ac es 3.4 
Pa) Saye ety) (5.4) 
Solving equation 3.4 for x(ty) by taking the inverse of P results in 


= Si ie es 


teas 
or) 
~~” 


W here: 
X(tg) is the nxl vector of the computed states at time to. 
ply (ty) is the nx! vector of the true states. 


ple. is the nxl vector of the error in the computed states. 


Defining 
Xa(to) = Poly, (tg) (3.6) 


and substituting this in equation 3.5 yields 


CJ 
| 
—— 


X(tp) = Xa(tg) + Pte, 


Or 


MUG Cala ple (Sro} 


The left hand side of equation 3.8 is actually the error in the computed state, 


e.. Substituting this notation in equation 3.8 gives : 
e.=Prle oe, 
C m 


Recalling what has been said about the dependency among the rows of P, it 
was stated that if the dependency was high the system was ill-conditioned and the 
outcome was large numbers for the rows of Pa Equation 3.9 shows that the error in 
the computed states is a linear combination of the measurement errors. For equal 
values of error in the measurements, Le., Clin ea eee the entries in the e, 
vector will be different according to the magnitudes of the entries in the corresponding 
rows of P. The larger the sum of magnitudes of the entries in a row, the larger will an 
error in Vq(tg) be reflected in x(to). 

To recapitulate , the upper-bound error criterion was defined as the reciprocal 


of the ratio of the error in the calculated state-variable to the observation error. 





l e. 
6 =e (3.10) 
obsi 
Sic nic 
cm 
Where: 


6 obs; iS the degree of observability of the ith state variable. 
€; 1 is the error in the ith measurement. 
e;, 1s the error in computing the ith state. 


Equation 3.9 can be rewritten as equation 3.11, where the ai; are the entries of P™!. 
C10 218m t 21 2"am Ts 7 Ela enm (Serr) 


C90 = 891 Ey mt 472 T+ T Bonfam 
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Substituting equation 3.11 in 3.10 vields 


e: 
ne (ani) 
7 Cea an Pes +G:_e 





ears, om 
obs 
(4 1€1m 
The main goal is to obtain a measure of the magnitude of the error in the 
computed states, for an equal error made while measuring the output. Therefore, it is 
reasonable to consider as equal all of the ey; [f one does not do so, misleading values 
for Onn for the state-variables mav be obtained, because the comparison between 
Bet Og e2: - O obs: Will be done under different conditions. 
Consider, then, that all the errors in the measurements are equal. Factoring 
out the e,, in the denominator of equation 3.12 and cancelling this value with the 


numerator, Vvields equation 3.13: 


I 


obs ys (a;5) 


§ (3.13) 


Since one is looking for the maximum possible value for the error in 
coniputing the ith state, the signs of the elements in the €4, vector must be such that 
an upper-bound for the computed error e;. will be obtained. Since the computed error 
1s a linear combination of the observation errors, the upper-bound for that error will be 


given by the sum of the absolute values of the row entries of P7!. This is expressed bv 


l 
= ried) 
obsi la ( 


Which is the final equation for evaluating the degree of observability for each state- 
variable using the upper-bound error method. 
2. Example 


Consider the single output linear time-invariant system described by: 


x(t)= Ax(t) + Bu(t) 
w(t) = Cxlt) 


17 


Whiere: 


me 7 B= *) c= fo 1] 


STEP ONE = Botner © smatiive 


[sss 
STEP TWO - Form on and normalize its rows, producing the matrix P: 


gt=fo 1 => p= [0 1 
PA -0.625 -0.78 | 


STEP THREE - Find P”!: 


p>! =| 1.249 -1.6 
L 60 


STEP FOUR - Find the degree of observability using equation 3.13: 


Q- er Q= : 4] 


| 


$88 = ee es 
Obst 171949) + ]-1.6| 

3- 
obs2 | seal) 


The upper-bound error criterion for x, says that the error on computing X(t) will be 
2.85 times bigger than the errors made while measuring y(tg) and its derivative if the 
assumption that both measurement errors are equal holds. The upper-bound error 
method criterion for X5 tells the user that the error on computing X>(tg) will be the 
same as the one made on the measurements. This result is verv satisfying because, 
looking at the C matrix, one can see that the state X 1s measured directly while x, 


influences the output measurement indirectly. 


D. STANDARD DEVIATION ERROR CRITERION 
1. Theoretical background 
Assume now that the errors in the output measurements that are being made 
are statistically independent random variables whose second order statistics are known. 
[t is Known that the variance of a linear combination of statistically independent 
random variables is the sum of the coefficients squared, multiplied bv the respective 


variances of each random variable [Ref. 11: p. 100]. Thus, taking the variance of the 


ith measurement and dividing it by the variance of the error in computing the ith state- 


variable gives 





im! 
Oi oc (set) 


If the stated assumption holds, i.e.. if the measurements of the output are 


statistically independent random variables, the following equation is obtained: 


Wate: s:| 
a — , —- (3.16) 
me Cacti oe at eye yee. Or, Nidlicarn 


eeneia 


Here again the assumption of equal variance is made because this is the only way of 
comparing the n degrees of observability under the same conditions. Factoring out the 


variances in equation 3.16 and cancelling with the numerator results in equation 3.17 


l 
a... = (ocda7) 
obs! 2 2 2 : 
Gy + E:5 +- agaoor aU 
| 
Ones vid (2 (3.18) 


Equation 3.18 is the final equation for evaluating the degree of observability 
for each state-variable using the variance as a measure of the error. 
2. Example 


Starting in step three of the previous exaniple: 


a ee inc 
1 0 


Find the degree of observability using equation 3.18: 
| 


6...) >= oe = 0.243 
obs] “(-4,249)- + (-1.6)" 
7) Peg ee l 
—— = 
obs2 Cis een 


The variance eriterion for evaluating the degree of observability shows that 
the error in computing Xx,(ty) has a variance 4.12 times bigger than the one 
encountered in the measurements. The variance of the error in evaluating X>(tg) 1s 


equal to the variance of the error in the measurements. 


FE. GRAPHICAL METHOD 
Another interpretation of ill conditioning 1s obtained by examining graphically a 


svstem) of two equations and two unknowns. Consider the following system: 


Q!x(to)=H(tg) (3.19) 
Wiens: 


Q! =! 1.01 0.99 
0.99 1.01 


Substituting the matrix Q Din equation >. lo iaekis 


L.01X (tg) + 9.99X>(ty) = ¥ } (ty) 


—~—_ 
U2 
tr 
a) 

~~ 


0.99X (ty) + 1.01X4(tg) = ¥a(ty) e2p 


Equations 3.20 and 3.21 can be seen as the equations of inom linesm ame 
abcissa X,(tg), the ordinate X(tg) and the x>(tg) intercepts at X)(tg)=9 being 
proportional to yj(tg) and y>(tg) respectively. The measurements of the output are 
¥, (tg) and ¥>{tg) and the main interest 1s to evaluate how much an error or disturbance 
made in the measurements (e,,) is “amplified” when the computation of the initial 
States takes place. In Figure 3.1 the lines corresponding to equations 5:20 and 3-2 iene 
presented. The lines intersect in the point X),(tg), Xo4(tg) which is the actual solution 
of the svstem of equations. Suppose that some errors are made while measuring the 
outputs vj(tg) and ¥5(tg). How is this fact going to change the position of the original 
lines? Hlow is the solution of the original svstem going to be perturbed? As in the 


development and discussion of Ablin’s method [Ref. 7], the errors in the measurements 


will be considered as equal so that the comparison between the different degrees of 
observability will be done under the same conditions. Consider, then. that the value 
measured for ¥ (to) exceeds its real value by an amount that is equal to €,,; In this 
case the line corresponding to equation 3.20 is displaced in the direction of positive 
Xr{tg) and X;(tg) to a position parallel to its original one. Conversely, if the measured 
value is less than its real one by an amount that is equal to € yy: the line corresponding 
to equation 3.20 is displaced to a position parallel to its original one but now in the 
direction of the negative X>(tg) and X, (tg). These two lines form an uncertainty region 


4 


around the line corresponding to equation 3.20. The same thing occurs for the line 
corresponding to equation 3.21 and a second uncertainty region is obtained. As can be 
seen in Figure 3.1 the intersection of these two regions forms a parallelogram shaped 
region that contains all the possible solutions of the svstem of equations if the 
measuremients can vary between ¥y(to)+ ey, and ¥j(to)-e,,, for the first measurement of 
the output and y5(tp)+e,, and v>(tg)-e,,, for the second one. The boundaries of the 
parallelogram: shaped regron contain the solution of the system of equations for the 
particular cases Where the errors in the measurements are maximum in magnitude. In 
Figure 3.1 eight different points are marked on the boundary regron. If the errors in the 
measurements are such that vj(tg) is measured at its minimum allowable bound. 
namely y)(tg)-e,, and ¥>(tg) is measured at its maximum allowable bound, point I 1s 
obtained. In the case of point 2 no error was assumed in the measurement of ¥)(to) 
while the other measurement was corrupted by +e¢,,. Point 3 depicts the case where 
both measurements were affected by the same positive error. In the case of point 4, 


¥ (tg) Was measured with no error but the error in ¥,(to) is +e. Point 5 iS. IN a Way, 


m’ 
symmetric to point I: y;(tg) ts measured corrupted by its maximum allowable error 
while y5(tg) is measured with its nunimum allowable bound. In the case of port 6, no 
error is assumed while measuring ¥)(tg), but the second measurement of the output 1s 
affected by an error equal to -e,,,. Point 7 is symmetric to point 3: both measurements 
of the output are affected by the same negative error -e,,. Finally. point § 1s the case 
where no error was present in the measurement of y>(tg) but y)(tg) was affected bv an 
error equal to -e,,. In all cases the soluuon of the system of equations 1s given by the 


projection of these points onto the axes Xj(tg) and X>( tg). 
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Graphical interpret 
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The discussion above gives enough information to answer the question of how 
the true solution of the svstem of equations is perturbed when errors in the 
measurements exist. As seen in Figure 3.1, the projections of points 1 and 5 onto the 
axes vield the most perturbed solutions. On the other hand, the projection of points 3 
and 7 vield solutions that are as close as possible to the true solution, assuming that, if 
anv errors are made in the measurements, they are equal to + €,,: Therefore, if the 
errors in the measurements are of opposite signs and equal magnitudes, the most 
perturbed solutions are obtained. In other words, if the errors in the measurements are 
as far as they can possibly be from one another the corresponding solution of the 
system of equations is also as far as it can be from the true solution. Conversely, if the 
errors in the measurements are equal both in magnitude and in sign, the least 
perturbed solution for this maximum error case 1s obtained. Thus, points 1 and 5 
correspond to the worst possible case when the errors corrupting the measurements are 
maximum while points 3 and 7 correspond to the best possible case when the 
maximum errors in the measurements are present. The projection of points number 1 
and 5 onto the axes X;(tg) and X5(ty) gives the maximum possible uncertainty region 
for the solution of the svstem of equations. As can be seen in Figure 3.1, these 
uncertainty regions are of the same order of magnitude for both states but much larger 
than the errors made in the measurements. 

Suppose that the Ql matrix produces lines that are perpendicular. As can be 
seen in Figure 3.2 the points 1, 5, 3 and 7 produce solutions that are nearly equally 
disturbed. Furthermore. the uncertainty regions in the solution of the system of 
equations are of the same order of magnitude as the errors in the measurements. 
Therefore, for the cases where the lines describing the system are orthogonal, the 
uncertainty regions in the solution of the system of equations are of the same order of 
magnitude as the errors made while measuring the output. 

Notice the similarities and differences between the cases depicted in Figures 3.1 
and 3.2. In both cases the maximum errors in the measurements are assumed to be 
equal. However, in the case of Figure 3.1 the errors in the measurements are 
enormously amplified when the computation of the states takes place, while in the case 
of Figure 3.2 the amplification is practically nonexistent. In both cases the uncertainty 
recions for both state variables are nearly equal yielding the conclusion that both state 


variables are equally observable. 
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Figure 3.2 Graphical interpretation of ill-conditioned systems-I[1. 
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What is behind the difference of behaviour in these two cases? One thing that 
rapidly comes to mind is the relative position of the lines describing the system: under 
study. In the case depicted in Figure 3.1, these lines have almost the same slope, 1-€., 
they are nearly parallel. In the case of Figure 3.2 the lines are orthogonal. 

In the cases studied so far both state variables had approximately the same 
degree of observability. What happens if the states have different degrees of 
observabilitv? Figure 3.3 depicts such a case. Suppose that the Ql matrix produces 
lines whose angles with the X,(to) axes are very small when compared with the angles 
between the lines and the X3(tp) axis. In this case the projection of the points | and 5 
onto the axes X)(ty) produces an uncertainty region that is much larger than the errors 
in the measurements. The projection of the same points onto the X>(tg) axes produces 
an uncertainty region that is much smaller than the one obtained for the other state 
variable and it is also of about the same order of magnitude as the error in the 
measurenient. Thus, the same error made in the measurements is enlarged a substantial 
amount when the computation of xj(tg) takes place while it remains unchanged when 
X>(tg) 1s computed. Therefore, according to the definition of relative observability one 
can sav that state X4(Ty) is relatively more observable than state X 1 (ty). Notice, 
however, that if by any chance the errors in the measurements are maximum but with 
equal sign ( points 3 and 7 ), the uncertainty region for Xj(tg) is very smal} while. for 
X4(tg), the uncertainty region doesn't differ much from the case of points I and 5. 
Therefore, the sensitivity of the solution to the variation of the errors ts also 
meaningful insofar as the relative observability is concerned: if the state is weakly 
observable it is also more sensitive to changes in the errors of the measurements. If the 
State 1s strongly observable it is not very sensitive to changes tn the errors of the 
measurements. 

The three cases described so far revealed two kinds of information: one imbedded 
in the relative position of the lines describing the system and the other imbedded in the 
position of these lines with respect to the coordinate axis. The first one determines the 
order of magnitude of the uncertainty regions of the computed state variables: lines 
nearly orthogonal imply uncertainty regions for the states of the same order of 
magnitude as the uncertainty regions for the measurements. Lines nearly parallel imply 
the possibilitv of at least one uncertainty region for the states much larger than the 


uncertainty region for the measurements. 
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Case of different degrees of observability. 
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The second information is the more important because it seems to be directly related to 
the evaluation of the relative degree of observabilitv for each state variable. In the 
cases depicted in Figures 3.1 and 3.3 the relative positions of the lines were nearly the 
same. However, their positions with respect to the coordinate axes were not: in the 
case depicted in Figure 3.1 the angles between the lines and both axes were about the 
same while in the case depicted in Figure 3.3 the angles between the lines and X 1 (tp) 
were verv small and the angles between X>(tg) Were almost 90 degrees. As previously 
discussed, in the case of Figure 3.1 both state variables were equally weakly observable. 
In the case of Figure 3.3 state variable X,(ty) Was Weakly observable while X4(tg) was 
relatively strongly observable. Thus, the angles between the lines and the coordinate 
axes are strongly related to the relative degree of observability of the state variables: if 
these angles are small ( lines nearly parallel to the axes ), the corresponding state 
variable is weakly observable. If the angles are nearly 90 degrees the corresponding 
state variable is strongly observable. 

Figure 3.4 shows the case depicted in Figure 3.1 deprived of several auxiary lines 
to enable a better understanding of the discussion to follow. Consider the lines 
perpendicular to those describing the system: the equation of the line orthogonal to 


the line described by equation 3.20 that passes through the origin is given by : 
1.01X4(te) = 0.99x j (tg) (322) 


The equation of the line orthogonal to the line described by equation 3.21 that 


passes through the origin is given by : 
1.01X (tq) = 9.99x9(tg) b 5928) 
Consider, now, the unit vectors u,; and us such that uy € X,(tg) and u, € Xx(ty). As 


can be seen in Figure 3.4 , the vectors v, and v5 belong respectively to the lines 


described by equations 3.22 and 3.23 and are given by: 
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v, = 1.0lu, +0.99u, (Se 
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Figure 3.4 Relation between the row entries of QI 
and the vectors v,; and V». 
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The components of the vectors Vv; and v5 are the entries in the first and second 
rows of ql or the first and second columns of Q. As can be seen in Figure 3.4 the 
angles between the vectors v) and v5 and the coordinate axes are nearly equal, i.e., 
$b, ~B>. This Figure depicts the case where both state variables have equal degrees of 
observability. Thus, equal angles between the vectors. whose components are the row 
entries of the Ql matrix, and the coordinate axes may imply that both state variables 
are equally observable. 

Figure 3.5 shows the case depicted in Figure 3.3 where state variable X7(tyy) was 
relatively more observable than state variable X 1 (tg). In this case both Vy and Ney Bbgs 
“closer” to X>(tp), 1.¢., 8) <O5 and y) < 7p. 

Analyzing and comparing the two cases just discussed, the importance of the 
information imbedded in the angles between the vectors, ( whose components are the 
row entries of the gl matrix ), and X;(tg) and X>(tg), can be understood: the smaller 
the angle between the normal vectors v) and v> and an axis of the state-space, the 
luigher the accuracy of computing the state corresponding to that axis, i.e, the higher 
the degree of observability for that state variable. The angles between the vectors and 
the coordinate axes are directly related to the magnitudes of their components along 
these axes: the larger the magnitude of the component along a particular axes, the 
smaller the angle between the vector and the axis. As previously discussed, the smaller 
the angle the higher the accuracy when computing that state. Recalling that the rows 
of Ql are the vectors relating X4(to) and X4(to) with ¥4(ty) and Va(ty). an easv way to 
check the relative observability of a state variable is to look for the largest principal 
axis component of each normal vector and then associate these components with the 
state Variables. This idea has been used as an algorithm for the evaluation of the 
relative degree of observabilitvy of a svstem and is one of the three algorithms presented 


in the Fortran program TOBS in Appendix A. 


F. A MORE GENERALIZED METHOD 

The graphical method previously described is quite helpful for low-order svstems 
because it requires few operations, is rapid and gives a good feel for the rank of the 
state variables as far as relative observability 1s concerned. When the svstem becomes 
larger and more complex the direction of minimum error, although perpendicular to the 
direction of maxinum error, may be difficult to find due to the increasing number of 


vectors in a higher order state-space. What if instead of checking the whole space 
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spanned by qi one finds a way of comparing one particular vector against the state- 
Space axes to obtain the relative degree of observability for each state variable ? 

The direction of minimum error convevs considerable information about the 
relative observabilitv of the svstem. Let W).W,...W, be the columns of Q normalized. 
In [Ref. 7: p. 9] the author forms an observability function called w which is a scalar 


that 1s given bv 


tas 
t— 
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where u represents a vector in the direction of minimum error. The length of this 


vector 1s constrained to be of unit length. i.e. 


ulus (3.27) 


Equation 3.26 can be written as follows : 


.28) 


ted 
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The quantity inside the brackets in equation 3.28 is actually non [Ref 12.3. 220% 


With this result in mind, this last equation can be written as 
w=u'(Q.Q, ju (3.29) 


The goal of this derivation ts to find a vector that gives the direction of 
maximum error, vielding a minimum for the observabilitv function y . To solve this 
problem, therefore, the gradient of wy with respect to u must be found. In this problem 
there 1§ a constraint that is given by equation 3.27. In general. any quadratic form of 
the tvpe x Px ( P positive definite ) is bounded below bv the product of the minimum 
eigenvalue of P and the length of x and bounded above by the product of the 
maxinium eigenvalue of P and the length of x [Ref. 13: p. 121]. For the particular case 
of equation 3.29, recalling that the length of u was assumed to be unity, this result can 


be written as 


(yd 
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min 1s the smallest eigenvalue of Q,Q, 


max iS the largest eigenvalue of Q.Q, 


The same result can be obtained using direct differentiation of equation 3.29 with 
respect to u [Ref. 7: p. 10]. In this procedure the Lagrangian multiplier formulation 1s 
used where the constraint is given by equation 3.27. The Lagrangian multiplier 


formulation is given by equation 3.31. 
oT qi Du -a(ulu-1y=0 (3.31 
F U (OO Ue (este) Ba 


The details on the differentiation of equation 3.31 can be found in [Ref. 14: p. 288]. 


The result of the differentiation is given by equation 3.32 
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Rearranging equation 3.32 and premultiplying both sides bv u° equation 3.33 is 


obtained. 
u'(Q,0, u=u hu= ay (3.33) 


This equation shows that the nunimum value of the observabilitv function w is equal 
to the smallest eigenvalue of ORME Therefore, the direction of maximum) error ( 
Where the observability function is minimum ) will be given by the eigenvector 
associated with the smallest eigenvalue. The algorithm for the evaluation of the 
relative observability for each state variable 1s, now, evident and is as follows: 

Lee om @ 7 Ons and multiply them together - call the product matrix QQT 

2. find the eigenvalues and associated eigenvectors of QQT 

3. find the smallest eigenvalue and associated eigenvector 


4. sort the components_of the eigenvector associated with the smallest eigenvalue 
In increasing order of magnitudes 


cay 


associate the sorted components with the respective state variables 


6. the smallest component is associated with the most observable state while the 
largest component 1s associated with the least observable state 


G. SUMMARY 

lise derer tie concept Of relative degree of observability was given as a 
measure of the accuracy involved in determining the initial states from the 
measurements of the output. The concept of an ill-conditioned matrix was presented as 
Well as the problem involved in inverting highly ill-conditioned matrices. Based on the 
definition of relative degree of observability given and the description of the problem of 
an ill-conditioned matrix. four methods for evaluating the relative degree of 
Observability were presented: the upper-bound error method, the standard deviation 
error method, the graphical method and the generalized graphical method. Examples 
were presented for the first two methods and algorithms were derived for the first. third 
and fourth methods. These algorithms are all imbedded in the same Fortran program 
(TOBS) presented in Appendix A. 


IV. INMPROVEMENT OF THE DEGREE OF OBSERVABILITY OF A 
WEAKLY OBSERVABLE SYSTEM 


A. INTRODUCTION 

In the last Chapter. methods for evaluating the degree of observabilitv for each 
state Were presented. A natural extension of the study of this problem is the search for 
a method of improving the observability of a particular weaklv observable state. In this 
Chapter methods for improvement will be presented. both for the single and multi- 


output cases. 


B. THEORETICAL BACKGROUND 

Before attempting to derive the method to increase the relative observability of a 
state variable, some important concepts of linear algebra need to be reviewed. 
Consider, then, the system represented by equation 2.1 where the square matrix Audis 
such that all its eigenvalues are distinct. In this case the matrix A has at least one 


eigenvector for each distinct eigenvalue that satisfies equation 4.1 
Ap=kp (4.1) 


Where pO is the eigenvector of the square matrix A associated with the eigenvalue 4 
(A is a scalar) [Ref. 10: p. G61]. Since the assumption is that the eigenvalues of A are 
all distinct, the n corresponding eigenvectors will all be linearly independent. Let 
P=(pj.--Pyl be the matrix of the eigenvectors of A. The matrix M=P! exists 
because the columins of P are linearly independent. Let m,! be thesrows of Mi ieam 
be shown that Alm =A.m or equivalently m;+A=).m,!) lee m;! are often called 
the left ( or row ) eigenvectors of A in contrast to the right ( or column ) eigenvectors 
of A [Ref. 10: p. 664]. Note, however, that the so called left eigenvectors of A are 
nothing but the regular eigenvectors of Allin [Ref. 10: p. 664] the author shows two 
relations that will be quite important in the derivation of the method described in the 
next section. These relations are described as being the spectral decomposition of A 


and are given bv 


A =v rpm! (2b 2) 
=PAM 
=\lAM 


res 2 Pe 
At=y fe pm; ! (4.3) 
=PA-M 
=\laem 


In general the spectral decomposition of A™ is given by 


a — ss Mo. f (4.4) 
= 
=\lAny 


A,}- The relations presented in the equations 4.2, 4.3 and 4.4 


will be verv helpful in understanding the method that is presented in the next section. 


Mocree = diag (Aj5h5,...,4} 


C. DERIVATION OF THE METHOD 
Consider. now, the transpose of the observability matrix defined in equation 2.3: 


C 
CA 
gl= CA? (4.5) 
ean 
Substituting the representation of the matrix A given by equations 4.2 and 4.3, yields 


equation 4.6: 


ere vl 
CM!AM 
g!=|cuM lato 7-6) 


cy lanl yy 
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Notice that in equation 4.6 the matrices C. Mo! and M are common to all the rows. 


Letting 
L= cul, (ans 


substituting equation 4.7 in equation 4.6 and post-factoring the matrix M, equation 4.8 


is obtained : 


L 
LA 
Qgl=|LA? | (4.8) 


LAD! 


The matrix L, is a row vector of dimension Ixn, because C is a 1xXn row vector ( in the 


I 


case of single output ) and M™* is a nxn matrix. Thus, L can be written as 


With this in mund and recalling that that A is a diagonal matrix whose elements are the 


eigenvalues of A, the matrix qi can be rewritten as follows : 


by a a L., 
Lip Lae ee 
[= Z 2 2 
US ie ae eee (4.10) 
| L 4 n-| Lex n-l is x n- 1 | 
el | eee n 


Looking carefullv at equation 4.10 it can be seen that the matrix Qi can be written as 


the product of three matrices : 


gl=vGm (4.11) 


Where V is the Vandermonde matrix of the eigenvalues of A: 


l EO | 
ae oe ee os 

Y= a >” ae _— (4.12) 
a ae yar 

G is a diagonal matrix whose elenients are the entries of L: 

L, Vi wens Q 

G=]0 Pre 0) (4.13) 
0 3:7 GA 


and M is the matrix of the left eigenvectors of A. 

With this procedure the matrix gl has been transformed into a product of three 
matrices. Since V and M depend only on the plant and G depends on both the plant 
and the C matrix. it seems to be G that needs to be studied and changed in order to 
improve the relative degree of observabilitvy of a particular state variable. 

It was mentioned in Chapter III that the degree of ill-condition of Qi 1S 
Important in deternuning the magnitude of the errors (e,) in computing the initial 
states once an error is introduced in measuring the output (¢,,). Therefore, it is 
valuable to have a quantitative measure of the degree of ill-condition of the matrix gl 
Peiiceito. p. 121] the author defines the condition number of a matnx W as the 


product of two matrix norms : 
cond(W) = || Wi] |W"! (4.14) 


If the condition number is small the matrix W is not ill-conditioned. If the matnx W 1s 
ill-conditioned the entries of its inverse will be large numbers, yielding a large condition 
number for W. However, this can also be true when the elements of W are small even 
in the absence of ill-conditioning. Multiplying the two norms has a normalization 


effect, so the condition number is large only for an ill conditioned system [Ref. 15: p. 


a 


[22]. Observe that the condition number cannot be smaller than 1, which corresponds 
to the conuition numiber ol thedentty. matrix [cers scent 

Equations 3.3 and 3.4 show that the vector of the measurements of the output 
consists of the sum of the actual values of the output [y,(ty)] with the error (em) 
introduced while measuring the output. In other words Valto) = yalto) > © ya iielaaats 


equivalent to 
em=Yalto) - Yaltg) (4.15) 


But ¥a(ty) = Q*x(t9) and ¥q(tg) = Q | x(to). By substituting these expressions in 


equation 4.15, equation 4.16 is obtained : 


em = Q'X(tg) - Q?Xaltg) (4.16) 
= Q![x(t9)-x,(to)] 


The error in the computed states is given by 
e.=(Q!y le, (4.17) 


Taking norms and recalling that, in general, |JAB]| S|{Aj] ||B||, yields for the particular 


case of equation 4.17 the result 


Heol S1MQ4) MI Nepal (4.18) 
Applying norms to the last part of equation 4.16 the following is obtained : 

He all S11Q7 II Heel (4.19) 
Notice that equations 4.18 and 4.19 constitute, respectively, upper and lower bounds 


on the magnitude of the error in computing the initial states. This fact is shown in 


equation 4.20 
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Applying the same reasoning to the equations relating Xa(tg) tO ¥,(ty), yields for 


Xq(tg) the upper and lower bounds shown in equation 4.21, 


| Vata) | 


Fol lsaltolll $ QTY] Ilyg(to)I (4.21) 


Which is equivalent to the result obtained in equation 4.22 


Te 
| eae 


<= = (4.22) 
Yalto ll — UXA(toll Il¥q(to)Il 


Qty" 








Multiplying this last inequality by the inequality 4.20 and recalling the definition of 
condition number given in equation 4.14, equation 4.23 is obtained for the bounds of 


the relative error in computing the initial states. 


| Hera! < ell <(cond oy Wem (4.23) 
(cond Q') |lva(to)ll  [IXg(to)l IIva(ty)l 


This last result shows that the relative error in the computed solution (e,) can be as 
large as the relative error in the measurements multiplied bv the condition number. Of 
course it can also be as small as the relative error in the measurements divided by the 
EOmaiion Mumber. SO, when the condition number 1s large, the error in the 
measurements gives little information about the accuracy of X(ty). Conversely, if the 
condition number is nearly unity. the relative error in the measurements 1s a good 
ieeasiire Of the relative error in comiputing the states (Ref. 15: p. 123], and the error in 
the measurements and the error in the calculated states are of the same order of 
miagnitude. The goal of this whole procedure is to improve the relative degree of 
observability of a weakly observable state. Therefore, if the condition number of ol! 1S 
mininuzed the relative error in the measurements will be of the same order of 
magnitude as the error in computing the states. The optinuum result will be obtained 
for a condition number equal to one. In this case the upper and lower bounds on e, 


Will coincide. The condition number of gl is given by 


Be 


cond(Q!)=|1Q "11 iQ’ YI (4.24) 
Recalling equation 4.11 where qi = VGM and substituting this in equation 4.24 vields 


cond(Q!) =||VGMI] |I(VGM)}] (4.25) 
=||VGMI| ||MtGlv-ty 


But the norm of the product of three matrices must not exceed the product of the three 


nornis, 1.e., 


cond(Q!) S]|VII [GI IMI Mo Gy ivy (4.26) 
<=(cond V\cond G){(cond M) 


Observe that V, Vel Moand M7! cannot be changed because they are totally 
dependent on the dvnamic characteristics of the plant. Therefore, only by minimizing 
the condition number of G will the condition number of Q! be mininuzed. The matrix 
G is a function of the entries of C (equations 4.7 and 4.13), which are the only parts of 
the svstem that the user can change. On the other hand the condition number of on 1S 
a function of the condition number of G (equation 4.26). Therefore. the relation that 
exists between the condition number of G and the condition number of Qi can be seen 
as a function whose independent variable and abcissa is the condition number of G 


and the ordinate ts the condition number of Ql, 
(cond Q!)=f{cond G) (sim) 


The goal of this procedure is to find the nunimum of the condition number of Qi by 
changing a specific entry in the C matrix. In other words, the goal is to find the 
minimum of f. In order to obtain an algorithm that accomplishes this, the function f 
must be studied in more detail. Both the abcissa and the ordinate are condition 
numbers of matrices which, by definition, can never be less than one. Only in verv 
special cases is the condition number of a matrix equal to one. The identity matrix and 
all its multiples are among those special cases. Notice, also that the first row of on 1S 
the C matrix in the single output case. Therefore, only when cy, 1s different from zero 


and all the other entries in the C matrix are Zero, will the possibility of obtaining a 


4) 


nunimum for the condition number of ql equal to one be likely to occur. In other 
words, in the majority of cases both abcissa and ordinate will be greater than one. The 
matrix G plays an important role in the development of the algorithm that minimizes 
the condition number of ql by changing a specific entry in the C matrix. In order to 
understand how the matrix G is used in this derivation a more detailed study of its 


entries will be done. The matrix G is a diagonal matrix whose elements are given bv 


j= 1.2.....n (4.28) 


where the Cj are Uiementniessosees@ Matrix ana the mj; are the entries of the M’! 
matrix. 

The norm used to compute the condition number may be anv matrix norm 
[Ref. 13: p.118]. In the present case the maximum column sum norm will be used to 


compute the condition number and will be defined as follows: 
WGI, = max) |g; j=1,2¢2,n (4.29) 


Notice that both subscripts of g are equal because the matrix G is diagonal. In this 
particular case the sum is equal to the only non-zero entry in the column. The norm is 


then given by equation 4.30: 


Gly = max|g; i= Geen (4.30) 


Suppose that the maximum element in G is gep where gee is given by equation 4.28. 


According to equation 4.30 the norm of G will then be equal to 


Gly = igeel ay 
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Since the G matrix is diagonal its minimum entry corresponds to the maximum entry 


in G7! which is also diagonal. Its norm is given bv equation 4.32 : 


| 
Gr !\]=——— j= 1.2... (4.32) 
min|g::| 
=a) 
Or 
I é 
| a) 
Where g,. is the smallest entry in G and is given by equation 4.28 . The condition 


number of G is, by definition, given by the product of the two norms given in 
equations 4.31 and 4.33. Substituting equation 4.28 for the values of g.. and gep in 


equations 4.31 and 4.35 vields, for the equation of the condition number of G 


n 


DdCyjMp 

ats 

= (cond G) (4.34) 
n 


Leys 
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The ratio given bv equation 4.34 is also the ratio between the largest and smattest 
eigenvalues of G (G is a diagonal matrix), which is another definition of condition 
number of a matrix [Ref. 16: p. 38]. The right-hand side of equation 4.34 is a constant. 
Let the inverse of that constant be equal to K. Equation 4.34 can then be rewritten as 


follows: 


n n 
KD icing = dicyjm (4.35) 
ae jal 


Since there is only one equation (equation 4.35 ) to provide the changes in the C 
matrix to improve the relative degree of observability of the system, a decision has to 
be made upon the entrv in the C matrix for which equation 4.35 will be solved. 
Assunung that at this stage one has already computed the relative degree of 
observabilitv of the system, one knows which is the least observable state. If that state 
Is state X,, equation 4.35 will be solved for C,, because this is the entry in the C matrix 


associated with the least observable state. Expanding equation 4.35 vields : 


K(cy);my¢ eas Ca) Oleg eeeerea) ilo | (4.36) 


Cppay ot ey mF... Fc) m,. 
Rearranging and collecting conimon terms. 
(Kniye m7 .)C] fear + (Knie gure ehprat te +(Km,¢ cheers O (4.37) 
Solving for c), equation 4.38 is obtained 


I 
ee Bic, B,#0 (4.38) 
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Where in general B= (Kme -m,,). At this point every tool is available to develop the 
algorithm that will provide the new entry in the C matrix corresponding to the 
nunimum condition number of Ql, 

The initialization of the algorithm is done by selecting an abcissa equal to the 
condition number of G of the initial system and an ordinate equal to the condition 
number of Ql of the same system. To define the interval where the niuninium is to be 
evaluated. another point has to be chosen. As previously mentioned, the condition 
numiber of a matrix 1s a quantity that is always greater or equal to one. With this in 
mind, the choice of the second point selected for the initialization of the algorithm 1s 
the one where the abcissa (condition number of G) is equal to one and the ordinate is 
the corresponding value for the condition number of Ql, The computation of the 
condition number of Qi given the condition number of G 1s a three step procedure 
that can be summarized as follows: 


STEP-1 : Given the condition number of G (abcissa) solve equation 4.38 for cy,. 
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STEP-2 2 Find the sie on by using equation 2.3 and transposing it. 

STEP-3 : Find the condition number of the new ql (new ordinate). 

A flow-chart of the algorithm with the initialization and munimization procedure 
1s presented in Figure 4.1. In the flow-chart the letter f corresponds to the definition in 
equation 4.27 and represents the three step procedure described above. X represents 
the condition number of Ge while i stepresentc tne concition mumncrol, Ql. After 
bisecting the interval a decision has to be made concerning where to search for the 
nunimum of f, te. if to the right or to the leit of N37 ihe inioemation civen sp sci 
derivative of f near X3,Y3 1s used to make this decision: if the derivative is positive the 
interval to be searched is to the left of X3. If the derivative is negative the interval to 
be searched is to the right of 3. Once the direction of search 1s Known ties perms 
linuting the search interval are renamed: the point with the largest ordinate is X1.Y1 
and the point with the smallest ordinate 1s X2.Y2. If the stopping criterion is met the 
algorithm stops. Otherwise, the bisecting procedure is executed again. This algorithm 
gives a local minimum. Other initializations have to be used to trv to find other 
nunimums. The whole algorithm containing the required nuninuzation of the condition 
number of gl, was implemented in a Fortran computer program called COBS and 1s 
presented in the Appendix B. 

The choice of the entrv of the C matrix associated with the least observable state 
as the one to be changed is made here. However, this issue 1s controversial because real 
cases may happen where the user cannot change the C matrix according to the least 
observable state. Furthermore, an independent change in an entrv of the C matrix may 
not be possible, 1.e., a change in one entrv of C may only be possible if another entrv is 
changed. Several combinations are possible. Therefore, according to a specific case, the 
user may have to try several allowable changes in the C matrix to determine the 
nunimum condition number of gl, The program that is in Appendix B has the default 
feature of changing the entry of the C matrix associated with the least observable state. 
However, any other entry in the C matrix can be changed in order to increase the 
degree of observability of a specific state variable. This is acconiplished by following 
the documentation contained in the listing of COBS. 

Once the condition number of on has been minimized and the new entries of the 
C matrix found, the system’s observability has to be checked. In equation 4.11 it was 
Semen tacit qi could be deconiposed mtcm tie prodier Vor Sthneomem ie cee liic 


Vandernionde and the matrix of the left eigenvectors are full rank matrices because by 
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assumption the plant has distinct eigenvalues. Therefore, if G is a full rank matrix gl 
will also be a full rank matrix and the system will be completely observable. For G to 


be full rank all its entries must be different from zero because G is a diagonal matrix. 


D. THE MULTI-OUTPUT CASE 
What has been said so far is only valid for the single output case. Since this is 
very seldom the case, the multi-output situation will be considered in this section. 
Consider, again, the linear time invariant system described by equations 2.1 and 
2.2 where C is a mXn matrix and m 21. The observability matrix Q is, for the general 


case of multi-output, given by 
g=(cl alc! calc! .. calyrf ch (4.39) 


where € is the rank of C (Ref. 17: p. 503]. This new observability matrix is nxp, where 
p=(n-€+1)m. Notice that if m=1 and €=1 (single output case), the matrix Q will be 
square. The fact that the matrix Q is no longer square, transforms a relatively simple 
problem into a much more elaborate one. Recalling equation 3.2 where the relation 
between the derivatives of the output and the states is given, if one attempts to solve 
that system of equations for x(tg) when Q is not a square matrix, the concept of 
generalized or pseudo-inverse must be introduced. Premultiplying equation 3.2 by Q 


vields 
(th) = QQ! y 4.40 
Qvglty) = QQ! x(tg) (4.40) 
Solving equation 4.40 for x(to). 


x(tg)=(QQ*) Oy gto) (4.41) 


The matrix (QQ!) !Q is called the pseudo-inverse of Q! [Ref. 17: p. 577]. The 
solution of equation 3.2 using the pseudo-inverse is an approximate solution in the 
least-squares sense [Ref. 17: p. 378]. Let Qq! = R. The problem is, now, finding how 
invertible the the R matrix is, i.e., how ill-conditioned the matrix R 1s. One possible 
procedure is the one used for the single output case. Thus, performing the spectral 
factorization of the matrix A in the transpose of equation 4.39 by substituting its 


representation given by equations 4.2. 4.3 and 4.4 yields equation 4.42: 


i, 
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initialize by selecting 


Xi=cond G 
Yl=cond Q! 
X2= 1 
Y2={U) 
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Figure 4.1) Algorithm to provide the minimum of the condition number of Q! 
lv changing a specr{rc entry mm C. 
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Ql =| cw lazh (4.42) 
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Notice, again, that as in the single output case, the matrices C, M7! and M are 
eormmniom tO ail the rows. Let 


L=cmM"! (033 
where the dimension of L is mxn and 
et jen Salita 
C5] C49 eerie oy) 
C = C3] C39 Siac C3 (4.44) 
Cm] Cm? ra emn 
and 
My 4 ae Min 
M59 Nena m> 
34 bien M37 (4.43) 
mL? ae Mn 





Substituting equation 4.43 in equation 4.42 and post-factoring the matrix M., equation 
4.46 is obtained: 
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L 
LA 
Ql=| La2 M 


i 


(4.46) 


Performing all the matrix multiplications in equation 4.46 except those involving the 


matrix M, the following representation for ql is obtained: 


Cl S12 eeesoe 
on O49 ssereaiats 
‘m1 a eaeent 
SrAy 01249 eevcee 
SoAy Oo ee 


eeeonee 


(4.47) 


where Gi are the entries of L and a the svstem eigenvalues, are the entries of the 


diagonal of A. The goal is to decompose the first matrix on the right hand side of 


equation 4.47 as the product of two matrices such that only one of them depends on 


the C matrix. The matrix gi can actually be written as the product of three matrices: 


Qi = vem 


where VY has a Vandermonde structure, 
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(4.48) 


Ay Ay seeces ive 
v= | Aj? a oe (4.49) 
A ue Ao soeees Pa f 


Each matrix formung VY is a mxm matrix of the form 


hel ee () 
ae =| 0 hel . ( (4.50) 
0 ee Ae 


Where 1, the V super-row index, is equal to 2,3,....n and j, the V super-column index is 
equal to 1,2,...,n. Therefore, V has in each of its super-columns (n-€+ 1) matrices with 
m rows and m columns. On the other hand, there are n submatrices in each super-row 
of V. Thus, the matrix ¥ has p rows and q columns where p=(n-€+1)m and q=mn. 


The matrix I is an mxm identitv matrix. The matrix G 1s a qxn block diagonal matrix 


Ly Dy ee 0 
0 i (ee 0 

peers... .22hceeee 52:25 (4.51) 
0 i ae ie 


Each matrix contained in G is actually an mx! columin vector corresponding to the 


columns of the matrix L given by equation 4.43. In general, each submatrix of @ is 


given by 
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L= | 3. i=1,2,....n (4.52) 


This decomposition of the qi matrix seems to be coniplicated and even difficult to use. 
Two examples will be presented to clarify the ideas. Suppose that a third order plant is 
given in which only one output is available, i.e., n= 3, m=1 and the rank of C is C= 1. 
The dimension of gl is pxn Where p=(n-€+ 1)m=3 and n=3; V is a pxq matrix where 
p= 3 and q=nxm=3; G is qxn where q=3 and p=3. In summary Ql, V and G are 3x3 
matrices. The matrix VY has in each of its columns n-£+1= 3 matrices with m= 1 row 
each. According to equation 4.49 the first row matrices aré identity matrices which 
have mi rows (in this case one row), 1.e., the first row elements are three ones. The 
second row elements are found by using equation 4.50 along with equation 4.49. Since 
one is looking for the second super-row of V, the index i is equal to 2 and the super- 
column index j is equal to 1.2 and 3. Therefore, the matrices that constitute the second 
super-row of ¥ have m=1I row and have entries 4;, 4 and 3. The third and last 


super-row matrices of VY are We oe and a 


l l l 
2 2 Z 
hy ho h3 


The matrix G is a qxn, 1.e., 3x3 diagonal block matrix whose entries are given by 


equation 4.52. Since m=1 the if are scalars and are equal to a Loi 


£3=)3- The matrix g will be equal to 
G =| 0 $1 0 (4.54) 
0 0 613 


As can be seen by equations 4.53 and 4.54, bv applying the spectral decomposition of 
gl for this generalized case, a result is obtained that agrees with the one that would be 


obtained if equations 4.12 and 4.13 were used. 


50 


Suppose, now that the system 1s third order with two outputs. [n other words, 
n=3 and m=2. The rank of C is the smaller of m and n, ie., €= 2. Thus gl Wil be a 
(ets | minim mialtix, 1:¢., qi will be a 4x3 matrix. This matrix will be decomposed in 
the product of three matrices V, G and M whose dimensions are as follows: V is a pxq 
matrix where p=(n-€+ 1)m=4 and q=mxn=6. The matrix G is qxn=6x3 while the 
matrix M is nxn=3x3. Each submatrix of V is mxm= 2x2. Therefore. VY has two super 
rows and three super colunins because that 1s the only way of obtaining a dimension of 
4x6 for the matrix VY. The submatrices in the first super-rows of YV are 2x2 identity 
matrices. Therefore, in the first super-row of V there are three (n) 2x2 (mxm) identity 
matrices. The second super-row is computed using equation 4.50 with 1 (the super-row 
index) equal te two and j (the super-columin index) running from one to three since 


there are three super columns in the VY matrix. 


Te ae 
LO el. Ol ie aml (4.55) 

a ee Ones 

hp 0 dy 0 2; 


The entries of matrix @ are given by equation 4.52. Since m= 2 the Z; are 2x1 column 


vectors. The matrix G will have three (n) of these columin vectors 1n its diagonal. 


ce 0 0 
63] 0 0 
G =| 0 612 0 (4.56) 
0 $29 0 
0 y $13 
0 0 £23 


If one multiplies equation 4.55 (V) by equation 4.56 (G) the result obtained will 
be equivalent to the one that would be obtained if equation 4.47 is used for the same 
values of n, mand €. 

In order to find the initial states the inversion of a matrix must be performed, as 
was mentioned at the beginning of this Section. The matrix R=QQ! is the one that 


has to be inverted as shown in equation 4.41. Thus, for the multiple-output case the 


condition number of R is the value that gives a measure of how observable the state 


variables are. Expanding R the following 1s obtained: 


R = Qql (4.57) 
= (VEM)! (vem) 
= M'g¢ly!l vem 


[In general, any quadratic form of the tvpe x! px (P positive definite) is bounded 
below by the product of the minimum eigenvalue of P and the length of x and bounded 
above by the product of the maximum eigenvalue of P and the length of x [Ref. 13: p. 
121]. For the particular case of equation 4.57 this result can be written as 
xii? sxTQQhx sa 


ee ae 
Neon max!! XII (4.58) 


where ae and he are the maximum and minimum eigenvalues of the square 
matrix R= QQ~. The theorem mentioned in [Ref. 13: p. [21] cam be used Seweralsiimes 
in equation 4.58 until upper and lower bounds are obtained only as functions of the 
maximum and munimum eigenvalues of the square matrices viv, cl¢ and MIM. 


Equation 4.59 presents what happens for the upper bound 


xi ultelivl penx (4.59) 
Xv vyix tM lig leyoxy 
T Tonicleen 
maxi? Vr G°G)||x° M* Mx] | 
(Ve (Gs GeV vila 


xOQ! x 


lA 


ma x! 


lA 


ma x 


lA 


hma Md x! maxX 


and an equivalent result can be obtained for the lower bound, 1.e., 


KQQ1K2 hnin(Z Wamrin( GG min(M | MIX (4.60) 


niin! brain 


Recalling the fact that in general the norm of the product is less or equal than the 
product of the norms, the three results presented tn equations 4.58, 4.59 and 4.60 imply 
that 


hmax(QQ') Shmax' V' Vhmax( Sh max! M iM) | (4.61) 


a 


and 
NtOQ!) 2 d min! VV) min(G | GY min( Mi) oe) 


In [Ref. 16: p. 38] the author mentions that the ratio between the largest and smallest 
eigenvalues 1s a measure of how ill-conditioned a matrix is. With this in mind and 


looking at the inequalities 4.61 and 4.62 it can be said that 
cond(QQ!) <feond(v! vycona(g! g)feond(M ! M)] (4.63) 


Observe that although the condition number of aq! is a function of viy. gig and 
ME M.. any improvement of the degree of observabilitv of the system can only be 
accomplished by changing the condition number of ole because this 1s the only part 
of the product of the three matrices that is not totally dependent on the plant. 
Assuming that the user is only allowed to alter the C matrix and recalling equation 
4.43, it can be seen that only by minimizing the condition number of gl¢ wall the 
condition number of aq! be minimized. The matrix gle is a nxn diagonal matrix 
Whose entries are actually the lengths of the column vectors of @ squared. Let the 


entries of cle be the gi 


=e 
ae ao 
_ - sas 
SP fai 
1=*l,2eecm (4.64) 
Se 2 
fnn < in 


Thus, the condition number of eae. as it is defined in [Ref. 16: p. 38], will be given by 
the ratio of the largest and the smallest entries in gle. The condition number of c'l¢ 


iS given by equation 4.65: 


m 
Y Sie? 
ae =mcanid gg) (4.65) 
562 
Cis 
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Where Sip is the largest element in the g'¢ matrix while ¢,, is the smallest element in 
the same matrix. Comparing equation 4.65 with equation 4.34 and recalling the fact 
that the numerator and denominator of equation 4.34 are respectively the largest and 
smallest entries in the matrix G (which is the equivalent to the matrix G for the multi- 
output case), it can be seen that starting from two different definitions of condition 
number, two equivalent equations for the same quantity are obtained. The entries of 


the matrix G are also sumimations and are given , in general, by 


n 
Gij == yeaa i=1,2,...m J=1,2,...,n (4.66) 
Dame 


Substituting equation 4.66 in equation 4.65 for the values of Cp and ¢., yields 


m on 

2 
» (dei pMne ) 
1=1 p=1 


= (cond Gg!) (4.67) 
mM 1) 


¥ (LeipMps)” 
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At this point, if one attempts to use the same method used in the single output case, 
namely solving equation 4.67 for the entry of C that is related to the state whose 
relative observability 1s to be improved, it can be seen that this task 1s far from being 
an easy one becauSe the equations that have to be solved are quadratics. Therefore, 
one of the possible solutions for this problem is the following: the first assumption is 
that the user, after studying the possibilities of the system under test, knows already 
which entries of C can be changed and how much can they be changed. Then only one 
of the entries of C would be changed in each run over the allowed interval. This 
procedure would have to be repeated for every entry in the.C matrix that can be 
changed. Meanwhile, using equation 4.67 the condition number of gle would be 
computed as well as the condition number of Qa For every change in the entry of C 
under test the relative degree of observability of the system would be computed. A plot 


of the condition number of gle against the condition number of Qaql for everv 
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possible allowable change in the C matrix would be an extremely valuable output. 
With this information and the knowledge of what states are more observable under 
different conditions, a choice of the best allowable combination of the entries of C 
could then be made such that the condition number of QQ! IS Nunimum or a specific 
state variable is more observable than the others. Notice the parallelism that exists 
between the derivations involving the multi-output case and the generalized graphical 
method developed in Chapter III to evaluate the relative degree of observabilitv of the 
State-variables. In both cases the matrix QQ! MAS IMvO medias =the ‘centers oe! the 
problemi: a search for the minimum eigenvector in the method of Chapter III and the 
search for the minimum condition number of QQ! inne smielnou presemied iets 
Chapter. The goals of both methods are somehow different although correlated 


because the second uses information from the first one. 


FE. CASE WHERE THERE IS A CONTROL INPUT 

In the beginning of Chapter III the concept of relative observability was defined 
as a measure of how accurately one can deternune the initial states from the 
measurements of the output at time t = tj, ty,...t,. The control input u(t) was set to 
zero because the effect of the input signal can always be determined if X(ty) can be 
found. Now, considering the existence of a control input u(t)#0, the concept of relative 
observability must be redefined as follows: measuring y(t) and u(t) at time t = t). 
t>..... t, one wants to be able to evaluate Which state components are most and least 
accurately determined from the measurements of both the input and output. It will be 
shown in this section that the fact that one is measuring the input signal as well as the 
output signal doesn’t alter in any way what has been shown to be valid for the case of 
mge— Or Consider. again, the LTI system described bv equations 2.1 and 2.2. 
Differentiating equation 2.2 n-1 times and substituting x for its value given in equation 
2.1, the foltowing 1s obtained: 

V(t = Grit) (4.68) 

v( A = CAX(ty) + CBu(ty) 
V(t = CAPX(tp) + St o) + CBu(ty) 
... CAN a aan “Bu(ty) + CABU ty) + CBU(ty) 


(m-a) 


hi cathy X(ty CAM “Bu(ty) + CAT? Bulty) +... + CBulty) 
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Notice that the arguments of the input and its derivatives is ty which means that the 
derivatives of the inputs are being evaluated at time t=tp. Equation 4.68 can be 


Written in the matrix form 


Penere: 
Y¥q(tg) 1s a nxI column vector of the output and its n-1 first derivatives, 
X(ty) 1s a nxl column vector of the state variables at time to 
$2 is a nx(n-1) lower triangular Toeplitz matrix that represents the impulse 


response of the svstem and is given by equation 4.70 


0 0 _ 0 
CB 0 Pe 0 

Q =| CAB CB . 0 (4.70) 
!cAnM-2p CAMB. CBl 


and ug(tg) is a Ix(n-1) column vector of the input and its n-2 first derivatives. 


Considering, now, the single output case, equation 4.69 is solved for x(to) 


X(tp) = (Q7Y by 4(to) - (Q1) 'Qug(ty) (4.71) 


Observe that the projection of an error made in measurements of y(t) and u(t) in 
computing the initial states is determined, again, by the degree of ill-conditioning of 
qi which is the only matrix that needs to be inverted in this process. The matrix 92 
causes no problem to because it can be determined very accurately (assuming that the 
identification of the plant has been done). Consider the example used in Chapter ITI 


to explain the graphical method where 
Ql =! 1.01 0.99 
0:99 Soll 
The matrix 42 was equal to zero because the control input was set to zero. Now the 


control input is different from zero which implies that Q is also different from Zero. 


Let 92 be equal to 
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G2= 10 
u 


Where alpha is any real number different from zero. Substituting the matrices on and 


© in equation 4.69 vields 


0.99%) (tg) + 1.O1X5(tg) = ¥4(tg) : Gugy(ty) (4.73) 


Equations 4.72 and 4.73 can be plotted in the X4(to) X>(tg) plane as was done in 
Chapter III Section E. However, just by inspection one can see that the slopes of the 
lines corresponding to equations 4.72 and 4.73 are the same as the slopes of the lines 
corresponding to equations 3.20 and 3.21 although a new member exists in equation 
4.73. This new member appears due to the existence of a control input and is only 
responsible for a change in the intercepts, 1.e.. for a change in the solution of the 
svstem of equations of the example in Chapter HI. Since the slopes are equal so are 
the angles between the lines. It was shown in Chapter II] that the angles between these 
lines are highly correlated with the degree of ill-conditioning of the g! matrix: the 
smaller the angles the higher the degree of ill-conditioning of the matrix Ql, Therefore. 
the fact that a control input exists doesn’t bring any changes in the basic ideas 
presented so far leading to the evaluation of the relative degree of observability of the 
State variables of a given system. The same thing happens for the improvement of the 
degree of observabilitvy of a weakly observable state. In the case of u(t)=0 the 
procedure was the factorization of the gl matrix as the product of three full rank 
miatrices such that only one was dependent on the C matrix. which was the only part 
of the system that the user is allowed to change. Then, the goal was to find a new 
entry in the C matrix such that the condition number of gi was munimum. Once a 
muninmium was found the relative observabilitv tests were performed as well as the 
evaluation of the new condition number of oO" InietinisSscase where a controlminput 
exists. the same procedure is followed and. furthermore, the results obtained will be 
equal to the ones that would be obtained if u(t)=0 because the observability matrix Q 


is independent of u(t). 


The case of multi-output with control input is. again. very similar to the mult- 
Output case with u(t)=0. Consider equation 4.69 where the matrix Ql iS NOsloucer a 
nxn matrix. In order to compute the initial states the direct inversion of gl cannot be 
made because on is not a square matrix . Premultiplving equation 4.69 by Q equation 
4.74 1s obtained 


Q¥ y(t) = QQ! x ty) + QQua(ty) (4.74) 

and solving it for x(tg) 
j= ixel.. T,-1 : 
(Chyler eee) Qu lty) (4.75) 


Observe that in this case the degree of ill-conditioning of aq! is the one that has to be 
studied instead of the degree of ill-conditioning of Ql. For the evaluation of the 
relative degree of observability the generalized graphical method is the most surtable 
one because the matrix involved in the derivation of the algorithm is the Qa! matrix. 
For the case of tmprovement of the relative observability of a particular state variable 
the problems that were presented in the Section D of this Chapter are also valid here 


because the square matrix QQ! is independent of the control input. 


FL. SUMMARY 

In this Chapter a method to tmprove the relative observability of a weakly 
observable state was developed for the single output regulator (u(t)=0) case. This 
method was based on the concepts of spectral factorization of a matrix with distinct 
ergenvalues discussed in [Ref. 10: p. 609]. The method has the goal of mininuzing the 
condition number of the ql matrix. An algorithm was derived and a Fortran program 
(COBS) presented in Appendix B was written. This program improves the relative 
degree of observability of the least observable state. Documentation 1s provided to 
enable the user to improve the relative degree of observability of a particular state- 
variable other than the least observable one. This method was then extended for the 
multi-output case where the concept of the pseudo-inverse of a matrix was used. The 
factorization of on into the product of three matrices Where only one was not entirely 
dependent on the plant was performed and examples were presented to demonstrate 


that this general case factorization contains the single output case. The method to 


CA 
CO 


improve the relative degree of observability of a particular state-variable was then 
presented and shown to be different from the single output case as well as much more 
elaborate. Finally, the case where a control input is present was studied both for the 
single and multi-output cases, leading to the conclusion that the fact that u(t)#0 


doesn't change in any way what has been said for the cases where u(t)= 0. 


=o 


V. ANALYSIS OF THE RESULTS 


A. INTRODUCTION 

[In the last two Chapters algorithms have been developed to test the relative 
degree of observability of a LTI system and to improve the degree of observability of a 
weakly observable state. As previously mentioned, these algorithms are implemented in 
two programs (TOBS and COBS) whose performance will be diseussed in this Chapter. 
The results obtained while running the programs will be evaluated and some 


observations will be presented. 


B. PERFORMANCE OF TOBS 

This program tests and ranks the relative degree of observability for each state 
variable using three out the four methods presented in Chapter II]. namelv the upper- 
bound error method, the graphical method and the generahzed graphical method. In 
the first phase of testing the program, the performance of the methods was investigated 
using examples already tested and used in [Ref. 7]. The Table below shows the values 
for the degrees of observability obtained. As ean be seen all the three methods agree on 
What are the least and most observable states. The matrices A and C of this system 


are: 


Acs] -12 ° 3 Cr ee al 


TABLE 1] 
TOBS OUTPUT FOR ABLIN S SECON DIOR VE aie aioe 


Upper-Bound Graphical Generalized 
Method Graphical 
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Table 1 shows that in all the cases X 1s the most observable state. This makes sense 
because the matrix C shows that X> Is the state variable that is being measured. Once 
TOBS was tuned and running as desired. several state-space representations of transfer 
functions with different C matrices were studied. The maximum order tested was order 
ten. As the order of the plant starts increasing the graphical method begins disagreeing 
with the other two methods in the ranking of the least observable state. However, the 
results of the upper-bound and generalized graphical methods agree with one another 
in almost all cases. This behavior can be understood if one recalls the discussion of the 
graphical method presented in Chapter III. This method relies on finding the largest 
component of the vector relating the states with the measurements. As the order of the 
svstem increases the surfaces representing the relations between the states and the 
measurements start being planes and hyperplanes where the graphical reasoning behind 
this simple procedure 1s no longer easv to apply. However, even for high order systems, 
the most observable state 1s almost always the same in all three methods. The big 
advantage of the graphical method is the small amount of computation required, 
enabling the user to obtain a quick determination of the relative observabilitv of the 
system in low order or reduced order models. Furthermore, due to its relatrvely easy 
graphical interpretation, it can be quite helpful in understanding the ideas that serve as 
the basis for this tvpe of method. Although agreeing with the generalized graphical 
method, the upper-bound error method suffers from the fact that it needs a matrix 
Inversion to obtain the required results. With this method, the evaluation of the 
relative degree of observability of the state variables cannot be obtained when the 
systeni is not completely observable because the inversion of the observability matrix 1s 
needed. In anv of the other methods this information can still be obtained even if one 
Or more states are not observable. The generalized method seems to be the best one 
because although it presents a little more difficult graphical interpretation, 1t doesn't 
need anv matrix inversion and, more importantly, it can be readily used to compute the 
relative degree of observability for each state variable in the multi-output case while 
the other two methods, as they are developed. need some changes before they can be 
used in evaluating the relative degree of observability for each state variable in the 
multi-output case. 

As was previously mentioned, several runs were made using state-space 
representations of transfer functions from order two to order ten for different 


combinations of the entries of the C matrix. Since a large amount of computer output 
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would probably not be of much help tn understanding the relative performance of the 


methods described in Chapter III, a summary of the results 1s given in Table 2. 


TABI SZ 
SUMMARY OF THE RESULTS OBTALV ED Us. Ge TIRE: DIP PEs 
| METHODS 


Order of 
the Plant 





4 93% oe 99% 
5 88% 87% 97% 
6 82% S270 oo 
i 75% Faro 98% 
S 71% 71% 96% 


In Table 2 column one contains the order of the plant. Columns two, three and 
four contain values that are a measure of the correlation between the graphical method 
(number one) and upper-bound error (number two) -Rj, -, the graphical and the 


generalized graphical method (number three) -R,3- and finally, between the upper- 
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bound error method and the generalized graphical method “Kee, (he wales uy Table 2 


were computed using equation 5.1 


n 
R: = —HL i. j=1,2,3 (5.1) 


where n is the order of the system (number of states), n, is the number of runs 
performed for the same value of n but for different entries of the C matrix and Bees 
the number of matches between the two methods insofar as the relative degree of 


observability is concerned. Notice that the product n n. represents the maximum 


r 
possible number of agreements. For instance, suppose that the order of the systen is 8 
and one performs 10 runs. The maxinium possible number of matches is 80. Flowever. 
if the two methods only agree in 60 (n= 60) out of the 80 possible, the agreement 
percentage will be 75° and this is the value that appears tn the corresponding columns 
of Table 2. Looking at this Table one can observe what was mentioned earlier tn this 
Chapter, namely the degradation of the graphical method for the evaluation of the 
relative degree of observabilitv as the order of the svstem increases. Actually, the 
disagreement between the graphical method and the other two is concentrated in the 
least observable states. In most of the cases the most observable states were the same 
within one level of ranking in all three methods. As can be seen in the fourth colunin 
Sietable 2. the agreenient percentage 1s slowly decreasing as the order of the system 
increases. This occurs due the wavy the values are computed. In fact, the values for the 
degrees of observability for the most observable states are, in most cases, of the same 
order of magnitude for high order systems. A slight change in one digit can produce a 
different ranking in the state-variables, which is interpreted by equation 5.1 as a 
imismatch, producing a decrease in the agreement percentage. However, as the order of 
magnitude of the degree of observability is the same. one can consider that the two 
niethods still agree as far as the relative degree of observability 1s concerned. Of course, 
this decision can only be made bv direct analysis of the computer output. To illustrate 
this characteristic behavior the following example is presented. Consider the sixth order 


Sesten) where 
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As can be seen in columns four and six, although theoretically the upper-bound error 
and the generalized graphical methods do not agree in the two most observable states. 
they are, in fact. equally observable because their degree of observability is of the same 
order of magnitude within the same method. This just shows that one must be careful 
in arriving at any kind of conclusion just based on the output of the computer 


program. 
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The third and probably more interesting part of this study was the test of models 
of real plants. The first model studied was the reduced order model of a SL-7 hull with 
three degrees of freedom. The state space representation was derived in [Ref. 18: pp. 


955.592) and.is shown below : 


Vv -) .64e-2 eae 0 V 
r = -1.35e-4 =5,87e-2 0 r 
‘le 0 | () SH 
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Two runs of this model were performed with two different C matrices. The results 
presented in Table 4. show the output of TOBS for C = [0 1 1] 


TABLE 4 
Mees OU lIrePPreR THE SL-7HUPL NIODEL -CASE I- 


eEcias ound Graphical Generalized 


Nfethod Graphical 





and Table 5 shows the output of TOBS forC = [0 0 1 }. 

As can be seen in Tables 4 and 5, in both cases the three tests gave v, the sway 
rate, as the least observable state, which agrees quite strongly with one’s intuition 
about the problem. Furthermore, in [Ref. 18] the author mentions that for some tvpe 
of studies the sway rate can be neglected in a reduced order model because its influence 


in the behaviour of the third order model 1s not verv strong. 


TABEES 
TOBS OLTPLTBOR THE SE HULE MOD Er-@ ra 


U pper-Bound Graphical Generalized 
Error 


Method Graphical 





The second model studied in this phase 1s a seventh order plant that is presented 
as example two of [Ref. 19], where the author is seeking the identification of the 
essential states to reduce the order of the model. With his test he finds that states x4 X65 
Xg¢ and Xz are the ones that contain more information on the donunant modes of the 
plant and are, therefore, the ones to keep in the reduced order model. Performing the 
three tests described in Chapter II] on his model, the results shown in Table 6 are 
obtained. The relative observability tests performed by TOBS show that the most 
observable states are X4. Xs. Xg and x7 (with special strength in the three last ones), 
Which agrees with the results obtained by the authors of [Ref. 19]. An interesting point 
to notice in this example is that there is a difference of one order of magnitude in the 
degrees of observability of X4, x; and x5 between the graphical and the generalized 
graphical method. If only the graphical method was available one could conclude that 
X4.X,, X> and X3 were all weakly observable. However, the upper-bound error method 
gives magnitudes for the degrees of observabilitv of the state variables that are of the 
same order of magnitude as the ones given by the generalized graphical method. 
Therefore, in this example only x3, and perhaps x, would be classified as weakly 


observable states. 


C. PERFORMANCE OF COBS 
The algorithm developed in Chapter IV for the single output case lead to COBS. 


a Fortran program. that improves the relative degree of observability of the least 
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observable state. Several runs were made of plants of order two up to order seven. For 
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a few cases of the seventh order plants the program took a long time to perform the 
minimization of the condition number of ot By relaxing the minimization stopping 
criterion the computing time decreased a bit but not dramatically. [lowever. the 
results that were obtained were verv satisfactory because the least observable state 
became the most or the second most observable state in all the cases with the 
corresponding minimization of the condition number of oe 

An interesting aspect to observe is what happens before and after the 
niinimization from the graphical point of view. Consider, as an example. the second 


order svstem presented in [Ref. 7] where initially 


A= 1-12 3 C= [0 Y 
Sed 


and the least observable state is X, as shown in Table |. Before the mininuzation the 


condition number of QI is 11.667 where 


gi =slo 41 
a 


After running COBS the new C matrix is C = [ 1.17 1 ]. Notice that the entry c,, 
was the one changed during minimization because this entry corresponds to the least 
observable state whose relative observability 1s to be improved. With this new C matrix 


the transpose of the observability matrix is 


Qh =fii7 1.0 

11.00 = -0.5 
whose condition number 1s now equal to 2.89. The new ranking of the state variables 
based on the relative degree of observability is shown in Table 7. Notice again the 
difference of one order of magnitude between the graphical and generalized graphical 
methods. The problem is now the inverse of what it was, 1e., now the graphical 


method agrees with the upper-bound error method in the order of magnitude of the 


degrees of observability. 


Pepe 7 
NEW RANKING OF THE STATE VARIABLES FOR ABLIN’'S EXAMPLE 


Upper-Bound Graphical Generalized 
Error Method Graphical 





The goal was reached, 1.e., the state that once was the least observable one is 
now the most observable and the condition number of Ql decreased. 
In order to obtain a graphical interpretation of what occurs in the system by 


changing the C matrix, the procedure used in Chapter III Section E can be used again. 
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Substituting the Q! matrix before the minimization in equation 3.19 the equations of 
two lines will be obtained. These two lines can be plotted in the X(t). Xo(ty) plane, 
assunung that ¥; and y5 are given measured values. For plotting purposes ¥y and y5 
are considered proportional to the two horizontal lines presented at the top of Figures 
3.1 and 5.2. Looking at Figure 5.1 it can be seen that the uncertainty region for X(t) 
is equal to the uncertamty region produced by the errors in the measurements. If one 
relates these results with those of Table | one can see that they agree because Rp isutine 
most observable state. Furthermore. two methods give for this state variable a value of 
| for the degree of observability which means that an error in the measurements is 
projected with the same order of magnitude when the computation of the initial states 
Pakeseplace. The state-variable X, 1s the least observable because its uncertainty region 
is larger than the one produced bv an error in the measurements. Substituting Ql after 
nuninuzation in equation 3.19 another different set of equations relating the states and 
the measurements of the output are obtained. These two lines are presented in Figure 
5.2 m the usual X (to). X>(ty) plane. The distances that correspond to Dead rate 
now the double of what they were in Figure 5.1 but their ratio is the same. This 
rescaling was done to obtain a cleaner plot of the two lines. Notice that now the 
uncertainty region for X;.,(tg) is of the same order of magnitude as the one produced 
by the measurements of the output while the uncertainty region for x5,,(ty) increased. 
The condition number of Ql decreased because the smaller angle between the lines tn 
Figure 5.1 ts about 33° and in Figure 5.2 is about 42°. 

In the runs of COBS with the different plants in about 40°%o of the cases the 
inttially least observable state became the second most observable, but with a value for 
the degree of observability of the same order of magnitude as the most observable 
state. Unfortunately, in the real world, the user may not be able to change the C 
matrix according to the least observable state. Generally there are constraints both for 
the entries of the C matrix that one is allowed to change as well as for the values that 
the entries can have. The program: COBS allows a change in the entry of the C matrix 
chosen by the user although the default is the the change in the entry of C assocrated 
with the least observable state. As previously mentioned, documentation tn the listing 
of the program ts provided to guide the user to set the program to change anv other 
entry in C. Despite the somewhat unrealistic philosophy upon which this algorithm was 
based. it provides a way of checking how the system and the states behave with respect 


to one another, which can be helpful to carry out further investigations. 
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D. SUMDIARY 

In this Chapter the performance of the methods used to evaluate the relative 
degree of observabilitv of a set of state variables was analyzed and some important 
results were obtained from the use of the program that performs these tests (TOBS). 
The analysis of the performance of the algorithm in improving the relative degree of 
observability of a weakly observable state was done and the graphical interpretation of 
a second order plant before and after the improvement was studied. Finally the 


inherent advantages and disadvantages of this algorithm were pointed out. 


de 


VineeINC LU SrhomN 


This thesis further develops the issue raised by Brown [Ref. 2] concerning the 
question of “how observable a linear system is”. This issue. which leads to the concept 
of relative observability. is addressed in four different ways in the thesis. Two of the 
methods used to test the relative observability of an LTI svstem are based on graphical 
interpretations of the relations that exist between the measurements of the output and 
the states. In the first graphical method the information imbedded in the angles 
between the vectors whose components are the row entries of the observability matrix 
Ode Te COOrdiInate axes (X(t). Xz(tg),--, Xp(tg)) Is USed 10 determine the melative 
Observability of the state variables. The smaller the angle between any vector and an 
axis of the state-space, the smaller the uncertainty region and the higher the degree of 
observability for that state-variable. Since the angles between the vectors and the axes 
are directly related to the magnitudes of their components along these axes, an easv 
algorithm to check the relative degree of observability of a state-variable was derived. 
constituting one of the three methods contained in the computer program TOBS. 
During the study of this method it was observed that a weaklv observable state was 
more sensitive to changes in the magnitudes of the errors in the measurements than a 
strongly observable state. The second graphical method is called the generalized 
graphical method and is developed from the idea of searching for the direction of 
maximum error vielding a minimum for the observabilitv function. This direction turns 
out to be given by the eigenvector associated with the smallest eigenvalue of the QQ! 
matrix. Therefore. the smallest component of the eigenvector associated with the 
smallest eigenvalue of aq! corresponds to the most observable state while the largest 
component corresponds to the least observable state. 

Two other analytical methods were presented and explained: the upper-bound 
error method and the standard deviation error method. Although these two methods 
seem to be conceptually different from the two graphical methods. a common thread 
was the fact that the degree of ill-conditioning of the on matrix plavs an important 
role in the theoretical development of both graphical and analvtical methods because a 
highly ill-conditioned matrix has a twofold meaning. On one hand it means that the 


rows of that particular matrix have a high degree of dependency (nearlv parallel) and 


on the other hand it means that the determinant of that matrix is a very small number 
(which causes large entries in the inverse of the matrix). The high degree of dependency 
among the rows of on is one of the basic ideas behind the development of the 
graphical method while the small number for the deternunant is one of the key ideas 
behind the understanding of both the upper-bound error and the standard deviation 
error methods. The graphical, upper-bound error and generalized graphical methods 
Were throughly tested in TOBS using different models. All three methods agreed in 
almiost all cases in ranking the relative observabilitvy of the states. However, the 
graphical method degraded as the order of the svstem increased. Despite this feature. 
the graphical method has a big advantage: it requires a verv small amount of 
coniputation, enabling the user to obtain a quick determination of the relative 
observability in low order and reduced order system models. Furthermore, due to its 
relatively easv graphical interpretation it can be verv helpful in understanding both the 
concept of an ill-conditioned matrix and the role that this characteristic of a specific 
matrix has in computing the relative observabilitv of a system. The upper bound error 
method 1s reliable but has two disadvantages: the requirement of matrix inversion for 
every evaluation of the degree of observability and the impossibility of using the 
method when the svstem 1s not completly observable. The generalized graphical 
method is the approach whose performance is more satisfactory because it can be used 
even if the system) is not completly observable, doesn’t degrade when the order of the 
svstem increases and doesn't need anv miatrix inversion for the required results to be 
obtained. Furthermore, it can be readily used to compute the relative degree of 
observability of a multi-output system while the graphical and upper-bound error 
methods require some changes before being readv to compute the relative degree of 
observability for each state-variable. 

In the second part of this thesis an algorithm to improve the relative degree of 
observabilitv for the single output regulator case was developed and tested in the 
computer program COBS. This algorithm was based on the concepts of spectral 
factorization of a matrix, allowing the decomposition of the ql matrix as the product 
of three matrices, two of them being fully dependent on the the plant and the other 
depending on both the plant and the C matrix. By altering the entry in the C matrix 
associated with the least observable state, the mininuzation of the condition number of 
on was achieved and the least observable state became the most or the second most 


observable state in all the cases that were tested. The single-output regulator case was 
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then generalized, first to a multi-output regulator problem and finally to a multi-output 
servo problem. The use of the concept of the matrix pseudo-inverse was needed for the 
theoretical development of the multi-output problem and the spectral factorization of 
gl was found to be possible although hard to accomplish. Two illustrative examples 
Were presented. An algorithm to improve the relative degree of observability of a 
specific state was also proposed, assuming the previous knowledge of what entries in 
the C matrix could be changed and how much could they be altered. Finally it was 
shown that the servo problem added nothing new to what had been accomplished for 
both the single output and the multi-output case because the kev matrices (Qi for the 
single output and aq! for the muulti-output) are always independent of the control 
Input. 

In conclusion, amongst the methods developed to test the relative degree of 
observabilitv of the system. it seems that the generalized graphical method is the most 
rehable one although the graphical method can give quick and easily computed 
information about the miost observable states. Since the upper-bound error method 
needs a matrix inversion, it may become computationally involved for high order 
systenis and therefore its use should be considered depending on the order of the 
system to be studied. 

By using the algorithms developed in this thesis. the designer can obtain much 
insight into a complex system with a relatively small amount of computations. Exactly 
how these ideas may be used will depend on the designer’s knowledge of the 
capabilities of the system as well as the specifications of the problem under study. 

There are two recommendations for further study. The first one is the 
iniplementation and test of the algorithm proposed for the improvement of the degree 
of observability of the system for the multi-output case. The second one 1s the study of 


the rssues addressed in this thesis for linear time-varving svstemis. 


APPEND 


TESTING THE DEGREE OF OBSERVABILITY OF AN LTISYSTENMI 


PROGRAM TOBS 


C 
CRAKKKAKAKAAKKRKAKARKA PURPOSE RARKAAKAKAAARKA AKA KAKAKARKARKA 


THIS PROGRAM EVALUATES THE DEGREE OF (@BSER Veco elec 
A SYSTEM, WHOSE MATRIX A AND @ ARE IN ADA Te eee ee 
THREE DIFFERENT METHODS: A UGRAPHI CARRIE tne Ur eee 
BOUND ERROR METHOD C AND THE GENERALIZED GRAPHICAL 
METHOD. THE Mak C ORDER FOR THE PLANT TSmeepere aie 
THE OUTPUT IS GIVEN € IN A DATACRILE WHER ow oom 
ARE PRESENTED IN DECREASING GRDER OF REBAT PP EVGroa. 
INPUT DATA IS FILE FrO4SFO0l)-— EGCICAL REECE Dei oec! 


OUTPUT DATA IS FILE FIOSrOO! = RGGI CALE SeGe. Deis se 
CAAKKAKAKK RK KKAR RK AK KK RK RAR RKA RK AK ARK ARK K ARK RR KARK KAKA 


MAANANANNANAN 
a tte A 


aM 


C*AKK* DECLARATION OF VARIABEE Ses = ~~ 
G 


IMPLICIT REAL*4 (A-H,0 
DIMENSION A(SO, Bee gt(50, 50) .Srii(50 3° TEST 50) DEOBs (50) 
DIMENSION IHARER 0} Cul oe 


DIMENSION WKAREA(130) ,ON(50,50) ,QQN(50, 50) 








COMPLEX W(50),Z(50,50) 
C 
C 
Cxx*** READ IN ORDER OF SYSTEM N, AND MATRICES A AND C *x*xx* 
C 
READ(4,*)N 
C WRITE(6,600)I 
READ(4,*)(C(1,1),I=1,N) 
DO 20 I=1,N 
C WRITE (6 700 
D(a. *) (ACL. a J=1,N) 
20 CONTINUE 
C 
C 
Cxx*A OUTPUT MATRICES A AND C ****x 
C 
WRITE(8,450 
WRITE(8,350)(C(1,1),I=1,N) 
WRITE(8,550 
DO 25 I=1,N 


WRITE(8,350)(A(I,J),J=1,N) 
25 CONTINUE 
C 


Cc 
CAAAAS COMPULE OC TRANSPOSE sae 
C 

DOwSsORI=15N 


ley =—C( led) 
30 CONT INUE 


Cc 
DO 60 I=2,N 
DO 50 JL=1,N 
Oni) — Or 
DO 40 JC=1,N 
OT(I , JLJ=OTCL , JL) +O teererenen al enna 
40 CONTINUE 
50 CONTINUE 
ae CONTINUE 
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/OTN(50,50) ,OTNI (50,50) , FLAG(50) 


e 
C*x*XXX NORMALIZE Q TRANSPOSE BY ROW **xxx* 


€ 
DO 80 I=1,N 
SUM=0. 
Do, Onge th 
SUM=SUM+QT(I,J)*OT(I,J) 
70 CONTINUE 


IF(SUM.EQ.0.)THEN 
QTN(I, JC)=OT (1 eo 
ON (IC, I )=Cniine 1c 
Po 782 yegntt, sey 
C)/SORT(SUM 
BHCC. ie Tae A ene 
75 CONTINUE 


ENDO IE 
80 CONTINUE 
C 


G 
Seta FIND BIGGEST ENTRY IN EACH COLUMN ***x*x% 
DO 100 J=1,N 


Ele —Or 
DO 90 I=1,N 
IF (ABS (OTN(I,J)).GT.BIG)THEN 
BIG=ABS(OTN(I,J)) 
END IF 
90 CONTINUE 
Be (d=). 


IF(BIG.EQ.0.)THEN 
EES 180053 
WRITE(8,800 
FLAG ( 

END IF 

END TE rg 

TEST(J)=STATE(J) 

100 CONTINUE 


C 
C 
ee COMPUTE THE RELATIVE DEGREE OF OBSERVABILITY ****x%x% 
EAnnXs UStNGelaeeUrrPER=BOUND ERROR. METHOD RARKKK 
C 

NSTOP=1. 

Doe 105s i=1., 


NSTOP= NSTOP*ELAG(I) 
105 CONTINUE 
IF(NSTOP. 8's .) THEN 
WRITE (6, 850 0} 
WRITE(8, 850 
END IF 
IF(NSTOP.EQ.1.)THEN 
CALL LINV2F(OTN,N,50,QTNI,0,WKAREA, IER) 
Domt20 1-1," 
DEN=0. 
Domo —20 
DEN=DEN+ABS(QTNI(I,J)) 
110 CONTINUE 
DGOBS(I)=1./DEN 
DGTEST(I)=DGOBS(I) 
120 CONTINUE 
END IF 


C 
CkARKR MULTIPLYING QN BY QTN ***** 
é 

DO 126e1=1.N 


DO124 JL=1,N 
QON(I,JL)=0. 


ay 


DO 122) ge=1, 
Q NCL, SL)= QON(I,JL)+ON(I,JC)*QTN(JC, JL) 
122 CONTIN 
124 CONTINUE 
126 CONTINUE 


Cs 
CxxxXKX FINDING THE SMALLEST EIGENVALUE AND CORRESPONDING RAKKK 
CRHAKKKAKKKKK EIGENVECTOR KAKKAKKKKK 


ena EIGRE (QQN, N,50,1,W,Z,50,WKAREA, IER) 
SMALL=1000 
DO 200 I=1,N 
IF (REAL(W(I)).LT.SMALL) THEN 
SMALL=REAL (W(1)) 
ISMALL=I 
END IF 
200 CONTINUE 
C 


C 
one PRINT IN DECREASING ORDER THE DEGREE OF OBSERVABILITY ***** 
CAA FOR EACH STATE VARIABLE AKRKKK 


WRITE(8, 880) 
DO 140 K=1,N 
BPRIN1=0. 
BPRIN2=0. 
BPRIN3=0. 
DO 130 I=1,N 
IF (TEST(L).GE.BPRIN1)THEN 
BPRIN1=TEST(I) 
IBIGI=I 


END IF 
IF (DGTEST(I).GT.BPRIN2.AND.NSTOP.EQ.1.)THEN 
BPRIN2=DGTEST(I) 
IBIG2=1 
END IF 
OTEST= (1.-ABS(REAL(Z(I,ISMALL)))) 
F(OTEST.GE.BPRIN3)THEN 
BPRIN3=OTEST 


IBIG3=I 
END IF 
30 CONTINUE 
c 
C 
CARKAKAAKKAKKARKKK OUTPUT *ARKAAAAKKKAAKKAKKR 
C 


PP Oe ETE (Oa THE 
WRITE(8 ,900)IBIG1, BPRIN1, 1BIG3,BPRIN3 
END WES 
TEC NSDOR: Te 1.)THEN 
ae Go; 910) IBIG1, BPRINI] IBIGZ(ERR 2 ease EER riIS 
EN 


TEST(IBIGL)=-1. 

DGTEST(IBIG2)=-1. 

Z2(IBIG3,ISMALL)= CMPLX(- Laat) 
140 CONTINUE 


e 
xAAK*X FORMATS KAKKKK 
C 
350 ‘FORMAT Gs. 13. ea Risa? 
370 PORMAT Qe bls. 77 4X Fee 
450 FORMAT (Uaioe 7 ee iis ENTRIES OF MATRIAvGCsARE om 
550 FORMAT(' ',//,' THE ENTRIES OF MATRIX A ARE :'! 
560 FORMAT(!3!)7/3""0 ITRANSP@sr tome 
570 FORMAT (ee TRANSPOSE NORMALIZED IS :') 
800 FORMAT(' ',//,'STATE X',12,1%,"IS NOP OSSHRVAEPLE es 
850 FORMAT(' ABLIN','S CRITERION CAN NOT BE EVALUATED, ',/,' BECAUSE 
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&Q TRANSPOSE IS SINGULAR! ,/) 

Soop Onwweee | RELATIVE DEGREE OF OBS. FOR GRAPHICAL ABLIN',' AND 
& GENERALIZED METHODS' ) 

900 eee ee eee 02 er oe 2's!) 2k,'K',12,F9.6) 

SOME OR tat meee IA 2aer eo, 20, XK! 12,2K,F9.6,2K,'X',12,2%,F9.6) 
END 


APPENDIX B 
INIPROVING THE DEG. OF OBS. OF A WEAKLY OBSERVABLE SYSTEM 


PROGRAM COBS 


c 
CARRKAKKKAKKKKKKKAAKKAKK PURPOSE KKRAKKAKKKAKKKAKAKARKARKKAKKAKKAKARK 


C THIS PROGRAM COMPUTES IE SRELATIVE DEGREE SOS ROR Beye ean: 

C OF THE STATE VaARIABERS OF A GIVEN SYSTEM By ieee e eo re eke 

C METHODS AND THEN VCOMECTIES THE NEW ENTRY SU ieee eee) 

C ORDER FOR THE DEGREE OF ILL~CONDITIONING One ItE Oboe vob led 1). 
C MATRIX TO BE MINIMUM: IN THES PROGRAM THE Se oro enoeee 

C CHANGED IS THE ONE CORRESPONDING TO THEVEBSS? Ob- ER SE cbs lore. 
C IF THE USER WANTS OR NEEDS 70 CHANGE OTRERGEIr yw sOR weer tae. 

C THAN THE ONE ASSOCTATED WITH THE LEAS TR@ESERVAE EES 02a tae 

C INTEGER VARIABLE LRVNEEDS TO°BE SET PORBioE ORDPReOhe enone 
C TO BE CHANGED. FOR INSTANCE, IF A4 NEEDS VTOSEE CHANGED isan ar 
C OF THE ENTRY ASSOCIATED WITH THE LEAST GBSERVABLE SI21E wns 

S 
e 
c 
C 
c 


USER SHOULD PUT LR=4 
KRRKKRKKKR ARK RRAKKKKAKKAKKKKKKKKKKAKA RRR AAR RR AR AK RAKARAAKKAKKRRA ARR 


kok OR OR OR FR OF 


AARKK YDECLARAl UOC ah Pae eomes a 
e 
IMPLICIT REAMCGeASH) OR, 


DIMENSION A(10, ate Re TEST(10) ,DGOBS(10) 
DIMENSION DGTEST(10),C(1,10) ,QTN(10,10$ ,QTNI(10,10) , FLAG(10) 
DIMENSION WKO(130) ,WK1(130) ,WR2(130) 
DIMENSION QN(10,10}, soon Pel tite ao 
DIMENSION V(10,10),G(10,10)},DA(10,10),DAR(10,10) ,QTI(10,10) 
DIMENSION AL(1,100),DARI(10,10) 
COMPLEX W(10),Z(10,10) 
C 
C 
Cxx***x READ IN ORDER OF SYSTEM N, AND MATRICES A AND C **x*x*xx 
C 
NT=0 
DEL=0.001 
READ(4,*)N 
READ(4,*)(C(1,1),I=1,N) 
DO 20 I=1,N 
READ (4,%*)(A(I,J),J=1,N) 
20 CONTINUE 
C 


c 
C*A*KX OUTPUT MATRICES A AND Meu came 


: 
WRITE(8,450 
WRITE(8,350)(C(1,1),I=1,N) 
WRITE(8,550 
DOn2s) T=1n 

WRITE(8,350)(A(I,J) , J=1,N) 
25 CONTINUE 
C 


C 
C*A**A COME UGE Os TRANSPO SH aes 
G 
CALESNULT (C7 Ay NO) 
c 
Cx*xxkk NORMALIZE Q TRANSPOSE BY ROW ***** 
G 
WRITE(6,700) 
700 FORMAT (! '' BEGINING COMPUTING THE DEGREES OF OBSERVABILITY' ) 


SO 


26 CONTINUE 


DO 80 I=1,N 
SUM=0. 
DO 70). J=1,N 
SUM=SUM+QT (I,J) *QT(I,J) 
70 CONTINUE 
IF(SUM. 80 Oe ) THEN 
OTN(I,JC)=0 or ee 
ON(JC,I)=OTN(I,JC 
DO 75 JC=1,N 
a oe ee eesti) 
ON(JC,I)=OTN(I,JC 
75 CONTINUE 


ENDOIr 
Boy CONTINUE 


ec 
—* FIND LARGEST ENTRY IN EACH COLUMN ****% 
DO 100 J=1,N 


BIG=0. 
DO 90 I=1,N 
IF(ABS(QTN(I,J)).GT.BIG)THEN 
BIG=ABS(QTN(I,J)) 
END IF 
90 CONTINUE 
aieG a) =1). 


IF (BIG.EQ.0.)THEN 
Hi ee 2 
WRITE(8,800)J 
FLAG(CE)=0. 

END IF 

STATE(J)=BIG 

TEST (J)=STATE(J) 

100 CONTINUE 
C 


ce 
C***X*X COMPUTE THE RELATIVE DEGREE OF OBSERVABILITY ****x%x* 
Cas UstiGetae  UPPER-SOUND ERROR METHOD AKKKKK 


C 

NSTOP=1. 

DO 105 I=1,N 
NSTOP= NSTOP*FLAG(I) 

105 CONTINUE 

Bra SHOE = EO een 
el et 
Pete a, clu 


END IF 
IF (NSTOP.EQ.1.)THEN 
CALL LINV2F(OTN,N,10,QTNI,2,WKO,IER) 
Dom loG raion 
DEN=0. 
DO 110 J=1,N 
DEN=DEN+ABS (QTNI(I,J)) 
110 CONTINUE 
DGOBS(I)=1./DEN 
DGIEST(1)=DGOBS(I) 
120 CONTINUE 
END IF 


c 
CARR S eee etnG ON BY QIN **<*~* 
C 


DO 126 I=1,N 
DO 124 JL=1,N 
OON(I,JL)=0. 
Homes 2esc=1N 
OON(I, JL)=QQN(I,JL)+ON(I,JC)*QTN(JC, JL) 


1 


ZZ CONTINUE 
124 CONTINUE 
ee CONTINUE 


ee FINDING THE SMALLEST EIGENVALUE AND CORRESPOND ING AAKKK 
CARKKKKKKAKKKA EIGENVECTOR KAKAKKKAAKK 
C 
CALL EIGRF(QON,N,10,1,W,Z,10,WKO, IER) 
SMALL=1000. 
DO 2Z0071- hi 
IF (REAL(W(1I)).LT.SMALL) THEN 
SMALE REAL(W(I)) 
ISMALL=I 
END IF 
200 CONTINUE 
€ 


E 
CxxxAX PRINT IN DECREASING ORDER THE DEGREE OF OBSERVABILITY ***x** 
CARKKK FOR EACH STATE VARIABLE aan 


WRITE(8, 880) 
DO 140 K=1,N 
BPRIN1=0. 
BPRIN2=0. 
BPRIN3=0. 
DO 130 I=1,N 
IF (TEST(1L).GE.BPRIN1) THEN 
BPRIN1=TEST (I) 
IBIGI=I 
END IF 
IF (DGTEST(I).GT.BPRIN2.AND.NSTOP.EQ.1.)THEN 
BPRIN2=DGTEST (I) 
nBIG2= 
END IF 
QTEST=(1.-ABS(REAL(Z(I, ISMALL) ))) 
IF (OTEST.GE.BPRIN3)THEN 
BPRIN3=QTEST 


IBIGS=2 
END IF 
1,510, CONTINUE 
S 
c 
CARKAKAAKAKKKAKKA QUTPUT AAKKAAKKKAKKKRAKKER 
c 
Dh CN Sere “ETE sly Welle 
WRITE(8, 360) 1BIGL, BPRIN], 1B1G3, BRhaie 
END TE 


IF(NSTOP. Ee .1.) THEN 
(WRITE 8,910)IBIG1,BPRIN1, IBIG2,BPRIN2,IBIG3,BPRIN3 
EN 


€ 
CRAKKKRKRKKAKKAKRKKKKAKKAKKAKKRKKKAKKARKAKKAKKAKAKKAKKAKKRAKKKAKKKRAKAKKKK 


CxARAX THE FOLLOWING THREE LINES OF CODE ARE THE ONES WHERE ****% 
Cxxx*x*X THE INTEGER VARIABLE LR IS ASSIGNED THE ORDER OF THE ***** 
Cxx*x* LEAST OBSERVABLE STATE. IF THE USER WANTS TO CHANGE ***x*x* 
Cxxxx* OTHER ENTRY OF C RATHER THAN THE DEFAULT ONE THESE 3 ****% 
CxxAKX LINES SHOULD BE COMMENTED AND THE 4TH ONE ACTIVATED ***xx 
CRAKKKKKKKKKKKAKKKAKAKKRKKKRKAKKKKAAKRAKAKKAKAAKKRAKKKRKAKRRAKAKRRKAAKKRRAKK 
C 

IF(K.EQ.N)THEN 

LR=IBIG3 
END 


Te 
C LR=ORDER OF IHE SIALTE (10 "BESGhAlcED 


TEST(IBIG1)=-1. 

DGTEST(IBIG2)=-1. 

Z(IBIG3,ISMALL)= CMPLX(- Lee) 
140 CONTINUE 


CALL CONDN(QT,N,N,CNQ) 
WRITE (8 ,610)CNQ 
WRITE(6,705) 
705 FORMAT(' ',' FINISHING THE DEGREES OF OBSERVABILITY' ) 
IF(NT.EQ.1)GO To 999 


C 
G2" *~78e2) 20th THE LEFT EIGENVECTORS OF A ***xxx 
C 

SREP EECRF(A,N,10,1,W,Z,10,WK1,I1ER) 


E 
Cx**x*X TRANSFORM W IN P WHERE P ARE THE REAL PARTS OF W **x*xx* 
C 


DO 1000 I=1,N 
DO 990 J=1,N 
P(1,J)=REAL(Z(I,J)) 
990 CONTINUE 
-1000 CONTINUE 


C 
Cx*x** FIND THE MINIMUM CONDITION NUMBER FOR Q TRANSPOSED ****x 
CxxAAAAAR AND CORRESPONDING NEW ENTRY IN THE C MATRIX **xxxxxxX 


e 
WRITE(6,710) 
710 FORMAT(' ',' MINIMIZING THE CONDITION NUMBER OF OQT') 
Y1=CNO 
GlreeMULTZ(C P,N-G,AL,LS,LG) 
LSOLD=LS 
LGOLD=LG 
CALL CONDN(G,N,N,CNG) 
Rileene 
c 
é 
CxxxRR COMPUTE THE CNQ FOR A CNG THAT EQUALS ONE ***** 
é 
PK=1. 
CALL CNEW(C,P,LR,LS,LG,PK,N,CN) 
C(1,LR)=CN 
CAEL MULT (C SE oe 
CALL CONDN(OT,N,N,CNO) 
Y2=CNO 
2 o 
C 
280 CONTINUE 
c 
E 
C*****k FIND THE ORDER OF THE LARGEST AND SMALLEST ENTRY IN G ****x 
é 
GAbeeMUETOCe P,N,G,AL,LS,LG) 
E 


Ce ey) ce 


Gee era xKX BISSEGT THE INTERVAL X1,X2 ***AAxRXX 
X3=(X1+X2)/2. 


CNG=X3 
E 
€ 
Cxx*** FIND NEW ORDINATE (CNQ) FOR THE PREVIOUS ABCISSA (CNG) ***** 
c 
PK=1./CNG 
CALL CNEW(C,P,LR,LSOLD,LGOLD,PK,N,CN) 
C( 1 9ER =cN 
CALL ee ,A,N, 1,97) 
CALL CONDN(OT,N,N,CNOQ) 
Y3=CNO 
E 


c 
Cxx*AX COMPUTE THE SLOPE NEAR X3 **x*** 


C 


C 

CxAAKK 
CARAKKAK 
CAAKKK 
CRARKK 


720 
Tze 
G 
G 


ADEL=X3+DEL 

CNG=ADEL 

PK=1./CNG 

CALL CNEW( G72, URyEseLD., Ler D erent sclE 
C(1,LR)=CN 

CALL SOUS Gna ,Aagtly yo) 

CALL CONDN(QT,N,N,CNQ) 


teepoie 
SLOPE=( svete 


RENAME THE NBWePOINTS SUCH THAT IN| THESINTERVAL 
TOWARDS THE DIRECTION OF DESCENDING Silose 
TRE POINT With THE LARGEST ORDINATE Sse ee 
AND THE ONE WITH THE SMALLEST ORDINATE Mise 2) 2 


IF(SLOPE.GE.O.)THEN 
IF(Y2.LT.Y1) THEN 
IF(Y3.LT.Y2)THEN 


mv 
ae 
Y2=¥3 
xe 
ELSE 
Y1=Y3 
an 
END IF 
IF(Y3.LT.Y1)THEN 
Y2=¥3 
ROSE 
ELSE 
v1 
Sen 
Y1=¥Y3 
a 
END IF 
END IF 
E 
IF(X2.LT.X1)THEN 
IF(Y¥3.LT.Y1)THEN 
Y2=Y3 
X2=K3 
ELSE 
Y2=Y1 
al 
Y1=¥3 
X1=K3 
END IF 
IF(Y3.LT.Y2)THEN 
i= 
K1=x2 
Y2=Y3 
X2=K3 
ELSE 
Y1=¥3 
=x 
END IF 
END IF 
END IF 
Thon os X1).GT.lE=5)Gompomece 
WRITE(6,720) 
FORMAT(' ',' OPTIMIZATION COMPLETED' ) 
WRITE(6,725) 
FORMAT(' ',' COMPUTING THE NEW DEGREES OF OBSERVABILITY') 


$4 


AARAAK 
KARKKK 
AKARKK 
AAARKAK 


. ’ ‘ 


C 


AAAN 


350 
370 
450 
500 
550 
560 
a0 
600 
610 


800 
850 


880 


900 
20 
999 


730 


e 


tas 50 


WEITE 


£300) pan wy 


WRITE (8, 600) 


Liss. 
NT= NTH 


IF(NT.EQ. 


1) GOTO 26 


KKKEKK FORMATS **x*x*x 


Fone ',F13. oe ee 

BOrtreet FILS. 7,/(4K,F13.7 

BOGE)" ,/7/>. THE ENDPRIES OF MATRIX C ARE :' 

FORMAT(' ',//,' THE NEW ENTRIES OF MATRIX C ARE :') 

FORMAT(' ',//,' THE ENTRIES OF MATRIX A ARE :') 

PORE: ',/je.) © TRANSPOSE IS :') 

FORME |) 7) eo TRAN GrOse NORMALIZED IS : 

BOREL!) eo tebe NEW DEGREES OF OBSERVABILITY ARE ;') 
este ie BUDE COMPL OOH NUMBER OF O TRANSPOSED IS EQUAL TO 

ie lee 

FORMAL IS) // wom tae i, | 1S) NOT OBSERVABLE! ,// ) 

FORMAT ABLIN','S CRITERION CAN NOE BE EVALUATED, ',/,' BECAUSE 
&® TRANSPOSE IS SINGULAR’ ,/) 

BORMAL( | “ST aeRELATIVE DEGREE OF OBS. FOR GRAPHICAL ABLIN',' AND 
& GENERALIZED METHODS ' ) 

See Po keer 2 ore Oa een 2 2K "Xk! 12 FS 26) 

FORMAT oe UR", 125255 eee 6 Les oy a FS, 6, Cee 2X, FS. 6) 
WeItE (6, 730) 

FORMAT(' ',!' PROGRAM COMPLETED’! ) 

SLOr 

END 


C 
CRAEKKKAKAAKKRAAAKKK SUBROUTINES XKAAAKKAAARARAKKKKKR 


C 
SeeececeoUBROULINEs tO COMPUTE THE COLUMN NORMS ***%*% 
Ee 


Gs 


THE PURO oMmMOoremaitos SUDROULINE IS 10 COnPUIE * 


C*THE MAXIMUM SUM OF THE COLUMNS OF A GIVEN MATRIX * 


& 


Cc 


SUBROUTINE NORM(DA,NL,NC,DBIG) 
REAL*4 DA(10,10) 


DBIG=-1. 
DO 3 J=1,NC 
TOC=0. 
note 
Poe=mOctaBs (DA( 1 .J)) 
CONTINUE 
LF (TOC.GT.DBIG)THEN 
DBIG=TOC 
END IF 
CONTINUE 
RETURN 


END 


c 
Cx**XX SUBROUTINE MULT1 (COMPUTES Q TRANSPOSED) **** 
oy 


Cx 


Pee rver@oneOr tHIS SUBROUTINE is Te COMPUTE a 


C*THE PRODUCT OF THE C AND A MATRICES TO PRODUCE QT” 
CARAAKKRAAAAKKKKKRAKAKAKRARKARARKAKRRKARAKAAKRKKARARKKA Kk 


SUBROUTINE MUET1(C,A,N,QT) 


CO 
ca 


REAL*4 C(1,10),A(10,10),QT(10,10) 
DO 30 Jada 
Teen = CCl wa) 
30 CONTINUE 


C 
Demeomi= 2a 
DO 50 JL=1,N 
ODT diye or 
DO) 40de se 
OT(1, JEV=0T (I, JL) +O f-eea erne 
40 CONTINUE 
50 CONTINUE 
60 CONTINUE 
RETURN 
END 


c 


Cc 
Ck*** SUBROUTINE MULT2 (COMPUTES THE MATRIX G) AND GIVES THE ORDER ***x 
CHAKAAKAR® OF THE LARGEST AND SMALLEST ENTRY IN THE G MATRIX ***%*xxx*x 


S 
CREAKKAKKAAKARKEKKK PURPOSE KAKAAKKRAKAAKKAAKARAKAKK 


C* THE PURPOSE OF THIS SUBROUTINE IS TO COMPUTE * 
C*THE PRODUCT OF THE C AND P MATRICES TO PRODUCE G * 
Cx IT ALSO PROVIDES THE ORDER OF THE SMALLEST AND * 
CX LARGEST ENTRIES IN THE G MATRIX AS WELL AS AL * 
C* WHICH IS A ROW VECTOR WHOSE ELEMENTS ARE THE * 
x ENTRIES OF G * 
CRRKKKRKKKARAKARRAAKKKRAKRRRKARKRRRARKARRARKRRARRRARARRKAKKA 
Cc 
SUBROUTINE ers Core LSenG) 
REAL*4 CCleaner 10,10) Aecimi0) Grnemon 
DO 1090 JL=1,N 
AL(1,JL)=0. 
DO 1080 JC=1,N 
AL(1,JL)=AL(1,JL)+C(1,JC)*P(JC,JL) 
1080 CONTINUE 
1090 CONTINUE 
GRL=-1.E+4 
POL=1.E+4 
DO 2020 J=1,N 
IF(AL(1,J).GT.GRL) THEN 
GRL=AL(1,J) 
G=J 


ENDS 

TF (An Gi eT. POL iHen 
tape 
Lo=J 


END IF 
2020 CONTINUE 
DO 1094 I=1,N 
DO 1092 J=1,N 
ere ted 
6a 


=f ae 
G(I,J)=0. 
END IF 
1092 CONTINUE 
1094 CONTINUE 
RETURN 
END 
Cc 
Cc 


cxxxx SUBROUTINE CONDN (COMPUTES THE CONDITION NUMBER OF AN ARRAY) **** 


SUBROUTINE CONDN (DAR,NL,NC,DCN) 
REAL*4 DAR(NL,NC) ,DARI(10,10) ,WK2(130) 
CALL NORM(DAR,NL,NC,BG) 

CALL LINV2F(DAR,NL,10, DARL 32 Wh2) ues 
CALL NORM(DARI,NL,NC,BGI) 


SO 


DCN=BG*BGI 
RETURN 
END 

e 


c 
G-“*~-~ SUBROUTINE TO COMPUTE THE NEW ENTRY IN THE C MATRIX ***x**% 
E 


SUBROUTINE CNEW(C,P,LR,LS,LG,PK,N,CN) 

REAL*4 C(1,10),P(10,10) 

ADD=0. 

DO 2040 I=1,N 
IF(I.NE.LR)THEN 

ADD=ADD+(PK*P(I,LG)-P(I,LS))*C(1,I1) 
END IF 
2040 CONTINUE 

CN=(-ADD) /(PK*P(LR,LG)-P(LR,LS)) 

RETURN 

END 
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